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Abstract 

Given a weighted graph, we introduce a partition of its vertex set such that 
the distance between any two clusters is bounded from below by a power of the 
minimum weight of both clusters. This partition is obtained by recursively merg¬ 
ing smaller clusters and cumulating their weights. For several classical random 
weighted graphs, we show that there exists a phase transition regarding the exis¬ 
tence of an infinite cluster. 

The motivation for introducing this partition arises from a connection with the 
contact process as it roughly describes the geometry of the sets where the process 
survives for a long time. We give a sufficient condition on a graph to ensure that 
the contact process has a non trivial phase transition in terms of the existence of 
an infinite cluster. As an application, we prove that the contact process admits 
a sub-critical phase on d-dimensional random geometric graphs and on random 
Delaunay triangulations. To the best of our knowledge, these are the first examples 
of graphs with unbounded degrees where the critical parameter is shown to be 
strictly positive. 
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1 Introduction 

The initial motivation of this work is the study of the contact process on an infinite graph 
with unbounded degrees. The contact process is a classical model of interacting particle 
system introduced by Harris in lU. It is commonly seen as a model for the spread of an 
infection inside a network. Roughly speaking, given a graph G with vertex set V, the 
contact process on G is a continuous time Markov process taking value in {0,1}^ (sites 
having value 1 at a given time are said to be infected) and with the following dynamics: 

• Each infected site heals at rate 1. 

• Each healthy site becomes infected at rate AN where A > 0 is the infection param¬ 
eter of the model and N is the number of infected neighbours. 

We give a rigorous definition of the contact process in Section]^ and refer the reader to 
the books of Liggett BZllS]] for a comprehensive survey on interacting particle systems, 
including the contact process. Durett's book [J^ also provides a nice survey on these 
models in the setting of random graphs. An important feature of the model is the 
existence of a critical infection rate Ac such that the process starting from a finite number 
of infected sites dies out almost surely when A < Ac but has a positive probability to 
survive for all times as soon as A > Ac. 

It is a general result that, on any infinite graph, there exists a super-critical phase 
i.e Ac < -|-oo (on a finite graph, the process necessarily dies out since it takes values in 
a finite space with zero being the unique absorbing state). This follows, for instance, 
from comparison with an oriented percolation process, see UZl. On the other hand, if 
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the graph has bounded degrees, then there also exists a non trivial sub-critical phase, 
i.e. Ac > 0. This can be seen by coupling the contact process with a continuous time 
branching random walk with reproduction rate A. Thus, the phase transition is non 
degenerated on any vertex-transitive graph such as and regular trees. The behaviour 
of the contact process on those graphs has been the topic of extensive studies in the last 
decades and is now relatively well understood. In particular, depending on the graph, 
there may exist a second critical value separating a strong and weak survival phase 
(see iri^lT3l[T4]| for such results on trees). Again, we refer the reader to ©IZilSl and the 
references therein for details. 

Comparatively, much less in known about the behavior of the contact process 
on more irregular graphs. Yet, in the last years, there has been renewed interest in 
considering these kind of graphs as they naturally appear as limits of finite random 
graphs such as Erdos-Renyi graphs, configuration models or preferential attachment 
graphs l|Til^|3ll4l[T0||Tl1l. 

However, without the boundedness assumption on the degree of G, the situation 
is much more complicated and the existence of a sub-critical phase is not guaranteed. 
For example, Pemantle HT^ proved that, on a Galton-Watson tree with reproduction 
law B such that, asymptotically, P{B > x} > exp(—for some j6 < 1, then Ac = 0. 
Thus, the degree distribution of a random tree can have moments of all orders and yet 
the contact process on it will still survive with positive probability even for arbitrarily 
small infection rates. This is a very different behavior from the one observed on regular 
trees with similar average degree and it indicates that the survival of the contact process 
depends on finer geometric aspects of the underlying graph than just its growth rate. In 
the case of Galton-Watson trees, we expect the critical value to be positive as soon as the 
reproduction law has exponential moments. This still remains to be proved, even for 
progeny distributions with arbitrarily light tails. In fact, to the best of our knowledge, 
there is no (non trivial) example of a graph with unbounded degrees for which it has 
been shown that Ac is non-zero. Worst, predictions of physicists hinting for non-zero 
values of A turned out to be wrong, see for instance HU. The main goal of this paper is 
to find a sufficient condition on a graph G for the contact process to have a non-trivial 
sub-critical phase and then give examples of classical graphs satisfying this condition. 

Let us quickly explain the difficulty met when studying the contact process on a 
graph with unbounded degrees. First, the comparison between the contact process and 
the branching random walk becomes useless since the later process always survives 
with positive probability. This follows from the fact that it can survive on finite star 
graphs with large enough degree. Therefore, we must find another way to control the 
influence of sites with large degree. Those sites should be seen as "sources" which, once 
infected, will generate many new infections. Some of those infections may, in turn, reach 
other sites with high degree. This can lead to an amplification effect preventing the 
process to ever die out depending on the repartition of these sources inside the graph. 
In this respect, it is not so hard to find conditions for the process to survive: one just has 
to find groups of vertices containing enough sites with very large degree to make the 
process super-critical. On the other hand, in order to ensure the death of the process, it 
is necessary to consider the global geometry of the graph. The following heuristic is 
meant to shed some light on this last statement and will motivate the introduction of 
a particular partition of the vertex set of the graph which we call cumulatively merged 
partition. This partition will, ultimately, become the main object of interest of this paper. 
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Figure 1: Illustration of the heuristic. In dotted lines are the maximal distances attained by infections 
started from each of the vertices a, b, c and d. Arrows represent infection fluxes and double arrows 
symbolize the grouping of sites, a group acting as a single site. In this example, a and c first merge 
together in (2), then the resulting cluster merges with site d in (3). Note that site b receives infections 
from the other sites but since it cannot reciprocate, it stays isolated during the merging procedure. 


Heuristic. It is well known that the contact process on a star graph of degree d ( i.e 
a vertex joined to d leaves) has a survival time of exponential order (say, to simplify 
exp(d)) when the infection parameter A is larger than some value Ac(d) > 0. Now, 
consider the contact process on an infinite graph G with unbounded degrees, and fix a 
very small infection rate A so that there are only very few sites in the graph where the 
contact process is locally super-critical (those with degree larger than say, dg). 

To see the influence of these vertices with anomalously large degree, imagine that 
we start the process with a single infected site a having degree da > do. In addition, 
suppose that in a neighbourhood of a, every vertex has degree smaller than do. Now 
run the process while forcing a to stay infected for a time of order exp (da) after which 
the whole star around site a recovers. By that time, roughly exp (da) infections will have 
been generated by the star around a. But, inside the neighbourhood of a, vertices have 
small degrees so the process is sub-critical and each infection emitted from a propagates 
only up to a distance with finite expectation and exponential tail. This tells us that the 
maximal distance reached by the infections generated from a should be roughly of order 

da- 

Now imagine that within distance smaller than da from a, there is some other vertex 
b with degree dh > do. Suppose also that di, is much smaller than the distance between 
a and b. The previous heuristic applied to b tells us that, in that case, the contact process 
started from site b has little chance to ever infect a. Thus, infections generated by 
a will propagate to b but the converse is false. This means that, while a is infected, 
infections regularly reach site b but this flux stops when a recovers, then b survives for 
an additional time exp{b) without reinfecting a. So, the whole process survives for a 
time of order exp{da) + exp{di,) ^ exp{da). 
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Consider now the case where there is a vertex c, again at distance less than da from a, 
but this time with degree dc also larger than the distance between a and c. In that case, 
infections generated by a can reach c and vice-versa. Consequently, when either site a 
or c recovers, the other site has a high probability to reinfect it before its own recovery 
This reinforcement effect means that, in order for the process to die out, both vertices 
a and c must recover almost simultaneously This will happen after a time of order 
exp {da) X exp(dc) = exp (da + dc)- Thus, for the purpose of studying the extinction 
time, we can see both vertices a and c acting like a single vertex of degree da + dc (see 
Figure Impart (1) and (2)). 

But now, our combined pair of vertices {a, c) will send infections to a larger distance 
da + di, and will possibly find other vertices to interact with (for example, vertex d in 
Figure]^. Iterating this procedure, we recursively group vertices together, with the 
condition that two groups merge whenever the sum of the degrees inside each group is 
larger than the distance between them. Assuming that this procedure is well defined 
and converges, the limiting partition should satisfy the condition that, for any two 
equivalence classes A and B, 

d(A, B) > min (r(A); r(B)} (1) 

where d is the graph distance and r(A) is the sum of the degrees of the vertices of A. 


In turns out that the limiting partition exists and does not depend on the order 
the merging procedure is performed. Its study is the purpose of Section where 
we rigorously define it for a general weighted graph. Then, we examine some of 
its properties and provide a description of its internal structure. To the best of our 
knowledge this cumulatively merged partition (CMP) defines a new model which 
appears to be quite rich while still remaining amenable to analysis. We think it is 
interesting in its own right and might prove worthy of further investigation. 

Next, in Section]^ we consider various types of classical random weighted graphs for 
which we study the question of percolation: is there an infinite cluster in the partition? 
We prove that, imder reasonable assumptions and similarly to classical percolation, 
there exists a phase transition. More precisely, we show that for i.i.d. weights, both sub- 
critical and super-critical phases exist on d-dimensional regular lattices (see Corollary 
|3.4[ Proposition |3.7| and Corollary 3^ for precise statements). We also prove similar 
results for geometric graphs and Delaunay triangulations weighted by their degrees 
(Proposition |3.12] ). Most of the proofs in this section rely on multiscale analysis as the 
recursive structure of the partition is particularly suited for this kind of arguments. 

In Section we return to our original question about contact process and we try 
to make rigorous the previous heuristic. The main result of this section, while not 
being completely satisfactory, gives a sufficient condition for the contact process to 
have a non-trivial phase transition by relating it with the existence of an infinite cluster 
for a particular CMP on the graph weighted by its degrees, see Theorem 4.11 As an 


application, combining this criterion with results from Section [33) . we obtain examples 
of graphs with unbounded degrees having a non trivial phase transition: 


Theorem 1.1. Let G be either a (supercritical) random geometric graph or a Delaunay trian¬ 
gulation constructed from a Poisson point process with Lebesgue intensity on Then, the 
critical infection parameter \cfor the contact process on G is strictly positive. 

Finally, in Section [^ we discuss some open questions and possible extensions related 
to the CMP and its cormection with the contact process. 
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2 Cumulative merging on a weighted graph 

2.1 Definition and general properties 

In the rest of this article, G = {V,E) will always denote a locally finite connected graph. 
In all cases of interest, the graph will be assumed infinite. We use the notation d(-, •) 
for the usual graph distance on G. The ball of radius I > 0 around a vertex x G V is 
denotedby B(xo,/) = {x G V,d(xo,x) < 1}. More generally, for A, B C V,we set 

d{A,B) = inf{d(x,i/) : (x,i/) G A x B}, 

B(A,/) = {z G V : d{z,A) < 1}, 
diam(A) = s\rp{d{x,y) : {x,y) G A^}. 

The graph G is equipped with a sequence of non-negative weights defined on the 
vertices: 

(r(x), X G V) G [0,oo)^. 

It is convenient to see r as a measure on V so the total weight of a set A C V is given by 

r(A) = Ya K^)- 

x^A 

We also fix 

1 < x < -|-oo 

to which we refer as the expansion exponent. For reasons that will become clear later, we 
call the quantity r[AY the influence radius of the set A. The triple (G, r, x) defines the 
parameters of our model. 

We are interested in partitions of the vertex set of the graph and need to introduce 
some additional notation. If C is a partition of V, we denote by the associated 
equivalence relation. In our setting, equivalence classes will often be referred to as 
clusters. For x G V, we denote by Cx the cluster of C containing x. Finally, If C and C 
are two partitions of V, we will write C <C' when C is coarser than C. The goal of this 
section is to study the partitions of V which satisfy the following property: 

Definition 2.1. A partition C of the vertex set V of G is said to he (r,x)—admissible if it is 
such that 


NC,C'^C, CfC' d(G,G0 > min{r(G),r(G0}'". (2) 

Remark. The most natural case corresponds to x = 1 when there is no space expansion. In 
this case, the definition of admissibility coincides with Condition^stated in the introduction. 
However, in later sections, we will need results valid for general x. Let us note that changing 
the expansion parameter from 1 to oc is not the same as merely changing every site weight to 
r(x)‘^. In fact, from the inequality (Er(x))'^ > I^r(x)“, we see that (r, oc)-admissibility implies 
{r‘^, l)-admissibility but the converse is false in general. 

Let us note that the trivial partition {V} is always admissible for any choice of (r, x). 
Given two partitions C and C, we can define 

Cnc' = {GnG' : G G C,C' G C'}. 
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More generally, given a family {C^)iei of partitions, we define Hfei via the equivalence 
relation: 


X 



c 

X ^ y for all i G 1. 


Suppose now that every partition is admissible. Let C ^ Che two distinct clusters of 
Pl/e/ C‘. By definition, there exist Zq G I such that C C C'o G C'o and C C C'o G with 
C'o Oo_ Hence 


d{C,C) > > min(r(C'°),r(C*°))“ > min(r(C),r(C))'^ 


which shows that plzei admissible. Thus, being admissible is a property 

which is stable under intersection of partitions. This leads us to the following natural 
definition: 


Definition 2.2. We call Cumulatively Merged Partition (CMP) of the graph G with respect 
to r and a the finest {r, a)-admissible partition. It is the intersection of all (r, a)-admissible 
partitions of the graph: 

^{G,r,cc)^ fl C. 

[r, a)-admissible 
partition C of G 

When there is no ambiguity, we will drop G, r,ccfrom the notation and simply write 

Remark. The CMP is monotone in r and x: for x < a' and for r <r' (for the canonical partial 
order), we have 

^(G,r,a) -<'Tf{G,r',x'). 

The formal definition of CMP is mathematically satisfying but it is not very useful 
in practice. We introduce another characterization which provides an explicit algorithm 
for constructing ^ by repeated merging of clusters (and justifies, incidentally, the name 
of this partition). 

Definition 2.3. Let x,y G V. The merging operator Mx,y is the function from the set of 
partitions ofV onto itself defined by 

A/r ^ ifCx^Cyandd{x,y)<mm{r{Cx),r{Cy)Y, 

1 C otherwise. 


i.e. it outputs the same partition except for clusters Cx and Cy which are merged together 
whenever they do not satisfy 

By definition, the merging operator always returns a coarser partition than its 
argument: 

C ^ Mx,y{C). ( 3 ) 

Moreover, it is easy to verify that this operator is monotone in the following sense: for 
any two partitions C, C and any x, y G V, we have 

C^C' Mx,y(C) -< Mx,y(C'). ( 4 ) 

Let us also point out that Mx,y restricted to the set of (r,x)-admissible partitions is 
the identity operator. We can now state the algorithm used to construct the CMP. It 
formalizes the procedure described in the introduction of the paper. 
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Cluster Merging Procedure (CMP). 

• Fix a sequence {xn,yn)neK of pairs of vertices ofV. 

• Start from the finest partition = {{x}, x ^ V} and define by induction 

• For every n, the partition is coarser than C”. Fhis allows to define the limiting 
partition 'if = lim t C” via the relation 

X ^ y X ^ yfor some n . 

Proposition 2.4. Assume that the sequence {xn,yn) satisfies 

{^n,yn} = {x,y} for infinitely many n, (5) 

for every x,y ^ V with x y. Fhen the partition ^ defined by the cluster merging procedure 
does not depend on the choice of the sequence {xn,yn) ond coincides with the CMP of Definition 

o 

Proof Let (C”) and (C”) be two sequences of partitions obtained by running the al¬ 
gorithm with respective sequences (x„,i/„) and (x„,y„), both satisfying property 
There exists an injection ^ : N —> N such that {xn,yn) = every n. 

Using the monotonicity properties Q and (|4]) of the merging operator, we check by 
induction that is coarser that C” hence ^ ^ and equality follows by symmetry 

The fact that the partition ^ obtained with the algorithm is (r, x)-admissible is 
straightforward since otherwise, we could find x, y G U belonging to different clusters 
and such that d{x,y) < mtn(r(^j;),r(^y))“; but then arid would necessarily have 
merged at some point of the procedure thanks to 

It remains to prove the minimality property Let C be an (r, x)-admissible partition 
and let (C”) be the sequence of partitions obtained by running the cluster merging 
procedure with the same sequence (x„,y„) but starting from the initial partition = C 
instead of the finest partition. Since, the merging operator does not modify an admissible 
partition and since ~< C^, the monotonicity of Mx,y yields 

C” ^ C” = C for all n. 

Hence, taking the limit n —^ oo, we get ^ C and so ^ is indeed the finest admissible 
partition. □ 

We start our study of the CMP by collecting some easy properties of this partition. 

Proposition 2.5. 

1. Any site x G U with r(x) < 1 is isolated in the CMP i.e. {x} G (the converse is false 

in general). 

2. For any C we have |C| < max(r(C),l). 

3. For any C G ^, we have the equivalence: r(C) = -|-oo |C| = -|-oo. 

4. Fhere is at most one infinite cluster. 
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Figure 2: [In this example G = N, a = 1 and radii not displayed are equal to 0] For a < 1, the CMP 
on this weighted graph is the finest partition (it contains only isolated vertices). For 1 < a < 1.5, 
there is a cluster with two sites and all the other sites are isolated. For a > 1.5, there is an infinite 
cluster composed of all the sites with non-zero radii. 


Proof. If a site x is such that r{x) < 1, then d{x,y) > 1 > r(x)'^ for any other site y, hence 
{x} will never merge during the CMP. This proves]^ Now, if a cluster contains more 
than one site, then according to[^, the radius of each of its site is at least 1, proving]^ The 
third statement follows from and the fact that all radii are finite by hypothesis. Finally, 
if there were two distinct infinite clusters, they would both have infinite influence radius 
so they would not satisfy Q. □ 

The proposition above states that the CMP remains unchanged if all the radii r(x) < 1 
are replaced by 0. On the contrary, changing, even slightly, the value of a single radius 
with r(x) > 1 can dramatically change the structure of This is illustrated in Figure]^ 
and it shows that local changes can propagate to infinity. 

Figure]^ also illustrates the fact that clusters of ^ are not necessarily connected sets. 
In tact, they can even have asymptotic zero density. The following proposition gives an 
upper bound on the diameter of the clusters with respect to their size and will prove 
useful in the next sections. 


Proposition 2.6. For any duster C G 

B{C,r{CY) isconneded 

and 


diam(C) < 

Proof. Define the function 

fix) = 


max 

r(C)'* 

2“-2 


hCiogzhc) 


2-,0j ifix = l, 

if a > 1. 


« J max ( if X = 1, 

ifx > 1. 


2^-2 


We first check that / satisfies the functional inequality 


( 6 ) 

(7) 


fia + b) > fia) -F fib) + min(fl, &)“ for all a,b >1. 
To see this, iixb > a > 1 and set z = b/a. We can write 


fia + b) — fia) — fib) — min(x, b)'^ = 


I (log2(l + z) +zlog(l + i) 




a 

2^-2 


( 8 ) 


2 ) if X = 1, 
if X > 1. 


In both cases, the functions appearing on the right hand side of the equality are non- 
decreasing in z and take value 0 at z = 1. Hence / satisfies 

Now, Q and Q are trivial when C is a singleton. In particular, all the clusters of the 
finest partition {{x},x G P} satisfy them. Since the CMP is obtained from repeated 
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merging operations starting from the trivial partition, it suffices to prove that @ and 
0 are stable by the merging operator. Assume that C\ and C 2 are distinct clusters for 
which ([^ and (0 hold and such that their merging is admissible: 

1 < d{Ci,C2) < min(r(Ci),r(C 2 ))“. (9) 

Then, clearly, B(Ci U C 2 , J'(Ci U C 2 Y) is connected. Moreover, using the triangle inequal¬ 
ity combined with (|^, we find that 

diam(Ci UC 2 ) < diam(Ci) + diam(C 2 ) +d(Ci,C 2 ) 

< /('•(Cl)) +/(''(C 2 )) + mm(r(Ci),r(C 2 ))“ 

</(r(Ci) + r(C2)) 

= /(r(CiUC2)). 


□ 


Remark. The bounds of the proposition are sharp. To see this when a. = 1, consider the sequence 
of weights An on {1,2,..., defined by induction by 


A^ = 11 and An+i = An 0... 0 A 


For each n, all the sites inside {1,2,..., n2” with radius 1 merge into a single cluster Cn 
with r{Cn) = 2" and diam(C„) = |2”. A similar construction also works for a > 1. 


2.2 Oriented graph structure on "rf 

The cluster merging procedure tells us that clusters of ^ are formed by aggregation of 
smaller clusters. This hints that there must be some structure hidden inside the clusters 
(because the merging operation cannot occur in any order). In this paper, we will not 
be concerned with this question even though we believe it could be of independent 
interest. Instead, we consider the relationship between distinct clusters of ^ and show 
that there is a hierarchical structure which provides a natural partial order over the set 
of clusters. 

Definition 2.7. We define the relation 1 —^ over the elements of^ by 

C^C' CYC' and d{C, C') < r{CY 

for any C, C' G YThen this holds, we say that C' descends from C. 

This relation is anti-symmetric: one cannot have C 1 —> and 1 —> C simultaneously 
because it would contradict the admissibility property of We interpret ('^, • 1 —> •) as 
an oriented graph on the set of clusters. The next proposition gathers some easy, yet 
important, properties concerning this graph. 

Proposition 2.8. The oriented graph (f€, • 1 —)• •) zs such that: 

1. The out-degree of any cluster C is smaller than \B{C,r{CY)\ — |C|. In particular, every 
cluster has finite out-degree except for the eventual infinite cluster. 
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2. If ^ has more than one cluster, then for every cluster C, 

C has out-degree 0 r{C) < 1 C is a singleton. 

3. IfC H-^ C', then [r{CY\ > [r{C'Y\ where [xj denotes the integer part of x. 

4. There is no infinite oriented path Ci C 2 ^ in (f€, 1 —(in particular there is no 
oriented circuit Q 1 —^ C 2 1 —> ... 1 —> Q). In addition, if C ^ ^ is such that |C| < 00 , 
then every oriented path started form C has length at most [r(C)“J. 

Proof The first statement follows from the fact that if C 1 — C', then rf(C, C') < r(C)“ and 
therefore C' Ci {B{C,r{CY)\C) f (Z). The second statement is straightforward recalling 
that G is assumed to be connected. For the third statement, notice that the graph distance 
d{-, •) only takes integer values. Consequently, if C 1 —^ C', then r{C'Y < d{C,C') < 
r{CY which implies < [j"(C)“J. For the last statement, if Ci i—> C 2 is 

a chain of cluster, then [r(C,)J is a strictly decreasing sequence in N U {+ 00 } so it is 
necessarily finite and its length is at most [j"(Ci)J (or [r(C 2 )J + 1 if Ci is the infinite 
cluster). □ 

According to the previous proposition, the graph , • 1 —)• •) does not contain any 
cycle. Thus, the relation [> given by 

C > there exists an oriented path from C to C' 

defines a partial order on ^ and Statement can be re-expressed in the form: 

Corollary 2 . 9 . Any totally ordered subset of ('^, l>) is isomorphic to one of the following 
ordinals: {1,.. . ,n}, N or N U { 00 }. 

Figure]^ gives some examples of possible oriented graph structures. It shows that, 
even though the out-degree of every vertex (except maybe one) is finite, in-degrees can 
be infinite. Let us remark that if there is one vertex with infinite in-degree, then there 
are infinitely many vertices with this property and we can find an infinite sub-graph 
H = {Ci,C2,...} CZ ^ such ths^t (T/ I— y C2j if 3.rici only if i ^ j. 


2.3 Stable sets and stabilisers 

In this section, we introduce the notion of stable sets. These are subsets of V such that 
the CMP on the inside and on the outside of the set are, in some sense, independent. 
Stable sets play a key role in understanding the structure of the CMP and in particular 
to determine whether or not there exists an infinite cluster. 

Let H be a subset of the vertex set of G. With a slight abuse of notation, we still 
denote by H the sub-graph of G induced by H (i.e. the graph with vertex set H and 
edge set obtained by keeping only the edges of G with both end vertices in H). We want 
to compare the CMP inside H, i.e. (H), with the trace over H of the CMP on the whole 
graph G which we denote by 

^{G)\h = {cnH: Ce ^(G), G n H ^ 0}. 

There is an easy inclusion: 
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Figure 3: Some examples of oriented graph structures on the set of clusters. In each example the 
underlying graph is Z and a = 1. Red vertices all have radius 1 and form an infinite cluster. Black 
vertices have radius 0. Blue vertices have radius written below them. In every example, all clusters 
are a composed of a single vertex except for the infinite cluster. 

First example: the oriented graph is finite. Second example: The graph is infinite but there is no 
infinite backtracking path. Third example: blue vertices form a infinite backtracking chain yet the 
in-degree of every cluster is finite. Fourth example: blue vertices form a infinite backtracking chain of 
vertices with infinite in-degrees [in this last example, arrows emanating from oo or pointing to black 
vertices are omitted]. 


Proposition 2.10. For every H C V, we have (H) -< ^(G)|^. 

Proof. Let {xn,yn)neiN be a sequence of pairs of distinct vertices of G satisfying We 
use it to construct '^(G). We can simultaneously construct ‘^(H) by considering only 
the indexes n such that (x„, yn) ^ arid the result follows from the monotonicity of 
the merging operator. □ 

Without additional assumptions, ^(G) can be strictly coarser than ^(H) since 
clusters growing outside of H can merge with clusters growing inside of H which, in 
turn, can yield additional merging inside H. 

Definition 2.11. We say that a subset H C V is stable (for the CMP in G) if 

B(G,r(G)“) C H. (10) 


Remark. 

• Being stable is a local property: we only need to look at the weights inside H to compute 

and check if it satisfies Thus, it does not depend on the value of the weights 
on G\H. 

• Since the CMP is defined as the finest (r, (x)-admissible partition, in order to show that a 
set H is stable, it suffices to find any (r, a)-admissible partition ofH satisfying (fTO]). 
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Figure 4: Examples of stable sets. Here G = Z and a = 1. The initial weights are displayed in grey 
above the line (black vertices are given weight 0). Clusters are each given a color and their zone of 
influence is materialized by an arrow. Some (but not all) stable sets are materialized by brackets (all 
the sites inside a bracket form a stable set). Note that the whole region displayed forms a stable set: 
the output of the CMP inside this region does not depend on the weights outside this region. 


The following proposition highlights the importance ot stable sets as its shows that 
they are the sets for which the CMP can be split into two separate partitions. 

Proposition 2.12. Let H C V be stable set. Then 

^(H) =^{G)\h and <^(G \ H) = <^(G)|g\h- 

Moreover, we have the decomposition 

^(G) = ^(H) U^(G\H). 

We point out that, even though H and G\H seem to play a symmetric role in the 
proposition, G\H is not necessarily stable when H is. A direct consequence is the 
following description of a stable set in terms of the CMP on G. 

Corollary 2.13. H C V is stable if and only if 


U Bix,r{^,iG)r) = H. 

xeH 


Proof of Proposition 2.12\ Fix a sequence {xn,yn) of pairs of vertices of G satisfying 
We start from the finest partition and simultaneously build three sequences of partitions 
(C^(G)) , and (C”(G\H)) using the merging procedure: 


C”+'{G) = M,.,,.(C”{G)); 

c«+t-(h) y« *= 

C” (H) otherwise; 

C"+i(G\H) = ifx„andy„ C G \ H, 

[ C”(G \ H) otherwise. 


( 11 ) 

( 12 ) 

(13) 


These sequences converge respectively towards 'if (G),'rf (H) and ‘^{G\H). Therefore, 
we just need to show that for each n, 

C^(G) = Cl (H) for every z C H, (14) 


and that similarly. 


Cl(G) = Cl(G\ H) for every z e G\H. 


(15) 
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Of course, we prove this by induction on n. Clearly, (14) and ( [T^ hold for n = 0 because 
we start from the finest partition. Assume that these equalities hold for n. If and 
yn are both in H or both in G \ H, then clearly the recurrence hypothesis still holds for 
n + 1. It remains to check that, if G H and i/„ G G \ H (the other case is symmetric), 
then no merging occurs in (3”"'“^(G). To see this, we compute 


d(Cl(G),C'l{G)) > d(CHG),G\H) 
= d(Cl{H),G\H) 
>d(V,,(H),G\H) 

>r{Cl(H)r 

= r{C’l(G)r. 


[rec. hypothesis] 
[H is stable] 

[rec. hypothesis] 


Thus d{xn,yn) > nairi(^"(C”^(G)),r(C"^(G)))“ which tells us that the clusters do not 
merge. □ 

Proposition 2.14. Let {Hi,i G I) be a family of stable subsets of V. Then 


IT = [^ Hi and ^ ~ 
iel iel 


are also stable. 

Proof Because each Hi is a disjoint union of clusters of ^(G), so are H and H. Fix 
C G '^(G). If G C H, there exists Iq such that C C Hi^ therefore B(C,r(C)‘^) C H/g C H 
hence H is stable. Similarly, if G C H, then B(C,r(C)‘^) C Hi for every i G I thus 
B{C,r{CY) C H and H is also stable. □ 

The following proposition provides a method for constructing large stable sets from 
smaller ones by "dilution": if a set is surrounded by stable subsets that are large enough, 
then the union of these sets is again stable. We will use this idea extensively in the next 
section to prove the existence of a sub-critical regime tor the CMP on random graphs. 

Proposition 2.15. Let W C W C V. Assume that 

W\W is a stable set 


and that 

B{W,r{WY) C W. 

Then W is a stable set. 


Proof. We simply observe that any cluster inside W has an influence radius bounded 
above by r(W)“ and therefore cannot reach outside of W. Since the subset W \ W is 
stable, its clusters carmot merge with clusters inside W and we find that 

U B(x,r(<^:,(W))") = U B{x,r{^,{W\W)r)U U B{x,r{^,{W)r) 

xew x€W\w xew 

C (W\ W) UB(W,r(W)“) 

C W 

which shows that W is stable. □ 
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The whole vertex set V is always stable. Since stable sets are stable under intersection, 
it is natural to consider the smallest stable set containing a given subset. 

Definition 2.16. For W C V, we call stabiliser ofW and denote by the smallest stable set 
containing W. It is the intersection of all stable sets containing W: 

sw^ n H. 

H stable 
WcH 

When considering stabilisers of single vertices, we use the notation Sx instead of 
Figure [5| shows some examples of stabilisers. We start by collecting some basic 
properties ofthese sets. 

Proposition 2.17. Stabilisers have the following properties: 

1. Sx = S^Jor any x G W 

2. D yJx(sW^x- 

3. S^^ = <Svv (ind ifWi c W 2 then <5^^ C 

4- = U;g7*Sw,- (ind Sn-^jWi C Cii^iSwi- 

5. For X G V, if Sx = ^x, then either r{x) < 1 and Sx = {x}, or Sx = V. 

Proof 1. and 2. follow from the fact that a stable set is a disjoint union of clusters of the 
CMP. Statement 3. is trivial and 4. follows from Proposition |2.14| and the minimality of 
stabilisers. The last statement is a consequence of the connectedness of G. □ 

Proposition 2.18. Let W G P. Define the growing sequence of subsets (S”) ofV by 

S0 = W and S”+^ = (J 

x€S« 


Fhen, we have 

= lim t S”. 

n 

As a consequence, ifC G then Sq is a connected set (although C itself need not be connected). 

Proof Let x G lim„ j" S”. Then, x G 5”° for no large enough hence B(^x,r(f^xY) £ 
5«o+i lim„ j" S”. This means that lim„ j" S” is stable. On the other hand, by a 
trivial induction argument. Corollary |2.13| gives that S” C for all n. Consequently 
lim„ t S” C Syq and equality follows by mini mali ty of stabilisers. The fact that Sq is 
cormected for C G ^ follows from Proposition 2.6| which ensures that is connected. 
Then S” remains connected for all n since we grow these sets by adding adjacent 
connected sets. □ 


A direct consequence of the previous proposition, stated in Corollary |2. 19 is that 
there is a nice interpretation of stabilisers of clusters in terms of the partial order [> 
defined in Section [Z^ It shows in particular that these stabilisers "pile over each other" 
as illustrated in Figure 
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Figure 5: Examples of stabilisers. The weighted graph is the same as in Figure]^ On the left, 
the stabilisers of each clusters inside the region are materialized by brackets. On the right, the 
corresponding oriented graph • i—>• •). Vertices and arrows pointing to black clusters of null radii 
are not represented. Notice that the leftmost cluster of size 1 is not a direct descendent of the cluster 
of size 8. 


Corollary 2.19. For any C we have 

5c = C u ( □ C') 

VcoC' / 

i.e. the stabiliser of a cluster C is exactly the subset composed of C together with all its 
descendants in the oriented graph (f^, ■ i—)• ■). Asa consequence, for any C, C' G we are in 
one of the three cases: 


Sc c Sqi, Sq o Sqi, Sc n Sc — 0- 

Moreover, we have the equivalence 

Sc = Sc ^ C = C. 


We can now state the main result of this section which will be instrumental in 
studying the existence of an infinite cluster in the CMP. The proof is straightforward 
but its usefulness promotes it to the rank of theorem. 

Theorem 2.20. [Criterion for the existence of an infinite cluster] Suppose that G is an 
infinite graph. For every x G V, the following statements are equivalent: 

1. = oo; 

2. \Sx\ = oo; 

3. <5;, = y. 


Equivalently, the partition ^ has no infinite cluster if and only if there exists an increasing 
sequence (Sn) of stable subsets such that lim Sn = V. 


C Sx, if Sx is finite then so is ^x- Reciprocally, suppose that '^x is finite, 
^x t> C} is a finite set (of finite clusters) according to Proposition 2.8 


Proof Since 
then {C G ^ : 

Consequently, Corollary |2.19| implies that Sx being the union of ^x and all these sets is 
also finite. Moreover, in this case we have Sx V since the graph is infinite. Finally, if 
^x is infinite, then Item of Proposition [Z5| asserts that its influence radius is infinite 
hence D B{^x,r{^xY) = V- □ 
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Remark. As we already noticed, changing the weight of a single vertex can have a macroscopic 
effect on the CMP. For example it can create or destroy the infinite cluster. Theorem 2.20 tells us 
that, on the other hand, there is still some locality in the CMP: clusters cannot "grow" from 
infinity. This essential feature of the model comes from the fact that we take the minimum of both 
radius in the definition of admissible partitions. Therefore, for two clusters to merge, both must 
reach to the other one. Most of the properties of the CMP described in this section would fail 
if, instead partitions satisfying 0 , we were to consider a definition of admissibility where the 
minimum of the radii of the clusters is replaced by either the maximum or the sum of the radii. 


We now give an algorithm that explores the graph starting from a given vertex Xq 
and reveals together with while never looking at the weight of any vertex 
outside of In particular, the algorithm stops after a finite number of steps if and 
only if the cluster containing Xq is finite. The fact that this algorithm works as intended 
is again a direct consequence of Propositions |2.12| and |2.18[ 

Stabiliser exploration algorithm. 

1. Start from the set of vertices = {xq} the partition = {{xq}}. 

2. At the n-th iteration of the algorithm, we have a finite set of vertices H" and a partition of 
this set C”. To go to the next step, we define 

ij B{z,r{C^y) andC^+^ 

(in practical implementations, the partition is obtained by running the cluster 
merging procedure starting from the partition C” U {{z},z G \ H”} instead of the 
finest partition on so that merges from previous iterations are not repeated at each 
step). 

3. If = H”, then the algorithm stops and outputs Sx^ = H” and ^ 5 ^.^ = C”. 
Otherwise, we iterate to step 2. 


Remark. Let us conclude our study of the general properties of the CMP by pointing out the 
fact that everything we established so far remains valid if we replace the expansion exponent by 
a general function i.e. if we consider partitions satisfying 

d{A, B) > i (min {r(A),r(B)}) 

for any pair of clusters A and B, where i is a non-decreasing function going toO atO and to +00 
at + 00 . The only notable difference is in Proposition \2.6\ where the upper bound for the diameter 
of a cluster is now given in term of a function f satisfying the functional equation @ with I in 
place of a. 

3 Phase transitions for cumulative merging on random 
weighted graphs 

In this section, we consider the CMP on several random weighted graphs. We investigate 
whether or not the partition ^ contains an infinite cluster. At first look, one might fear 
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that this will always be the case due to the amplification phenomenon resulting from the 
additive nature of cluster merging. Or, on the contrary, the cumulative effect might be 
quite weak and percolation by cumulative merging could be very similar to classical site 
percolation. It turns out that both worries are unfounded and that, for a wide variety of 
random weighted graphs, there is a non-trivial phase transition differing from that of 
classical site percolation. 

In this paper, we will consider the following three general families of graphs: 

Model 1 (Bernoulli CMP). The underlying graph G is a deterministic infinite graph (e.g. 

, a tree ...) and the weights {r{x),x G V) are independent identically distributed Bernoulli 
random variables with parameter p G [0,1]. We denote by the law of (G, r). 

Model 2 (Continuum CMP). The underlying graph G is a deterministic infinite graph and 
the weights {r{x), x ^ V) are independent identically distributed random variables with law 
AZ, where A > 0 and Z is a random variable taking value in [0, oo). We denote by the law 
of{G,r). 

Model 3 (Degree-weighted CMP). The underlying graph G is a random infinite graph (e.g. 
a Galton-Watson tree, a random Delaunay triangulation, a random planar map,...) and the 
weights are defined by r{x) = deg(x)l|jgg('j.)>^}, with A > 0. We denote by P^ the law of 

(G.r). 

The first two models are the counterparts in the context of cumulative merging of 
classical (site) percolation and boolean models. The third model may seem artificial at 
first. However, as we already explained in the introduction, it appears naturally in the 
cormection between cumulative merging and the contact process. We will investigate 
this relationship in Section 

For all these models, we will simply write P for the law of the weighted graph when 
the indices are clear from the context. Each model has a free parameter (p for Bernoulli, 
A for continuum, and A for degree-weighted CMP) so we ask, the expansion exponent a 
being fixed, whether of not ^ contains an infinite cluster depending on the value of this 
parameter. By monotonicity of the CMP with respect to a. and the weight sequence r, the 
probability of having an infinite cluster is monotone in both a and the free parameter of 
the model. 

Definition 3.1. Tor Bernoulli GMP we define 

Pc{oi.) =tnf|pG [0;1] : {‘^{G,r, a) has an infinite cluster} > o| G [0;1]. 

Similarly, for continuous GMP we define 

Xc{(x) = inf |a > 0 : P^'^ {^(G,r, a) has an infinite cluster} > o| G [0; +oo], 
and for degree biased GMP, we set 

Ac(x) = sup |a > 1 : P^ {'Z’(G,r,x) has an infinite cluster} > o| G [[l;+oo]] • 

Under fairly general assumptions on G, it is easy to check that the existence of an 
infinite cluster is an event of either null or full probability. 
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Proposition 3.2. Suppose that G is a vertex transitive graph, then for Bernoulli CMP (model 
1) or continuum CMP (model 2), we have 


P has an infinite cluster} G {0,1}. 


Proof Similarly to i.i.d. percolation on transitive graphs, the result follows from ergod- 
icity since the family of weights is invariant by translations and so is the (measurable) 
event of having an rnfinife clusfer. □ 


In the case of Model 3, one needs, of course, to make some assumptions on the 
random graph G in order to get a 0-1 law. However, for a large class of graphs, the 
existence of an infinite cluster is still a trivial event thank again to general ergodicity 
properties. This is in particular the case for the random geometric graphs and the 
Delaunay triangulations considered in Section 3.3 as well as for Galton-Watson trees or 
unimodular random graphs. 


3.1 Phase transition on Z for Bernoulli CMP 

By trivial coupling, it is clear that for any graph the critical parameter pc for Bernoulli 
CMP is smaller or equal fo the critical parameter psUe for classical i.i.d. sife percolafion. 
We now prove fhat these parameter differ in general. The following result shows that, 
even in dimension 1, there exists an infinite cluster in ^ for Bernoulli CMP when p is 
close enough to 1. This contrasts with the case of site percolation where psUeC^) = 1- 

Proposition 3.3. Consider Bernoulli CMP on G = N. For any ci>l,we have Pc{oi) < 1. 

Any infinite connected graph contains N as a sub-graph. Thus, by coupling, the 
proposition above shows that there always exists a super-critical phase for Bernoulli 
CMP on any infinite graph. Let us also remark that, if Z is a non-negative random 
variable which is not identically zero, then, for any 0 < p < 1 there exists, 0<A(p)<oo 
such that A(p)Z stochastically dominates a Bernoulli random variable with parameter 
p. Putting all this together, we get 

Corollary 3.4. Assume that G is an infinite connected graph and let a >\. 

• Pc{oi) < 1 for Bernoulli CMP. 

• Ac (a) < OQ for Continuum CMP as soon as Z is not identically 0. 

Proof of proposition |3.3[ Since the CMP is monotone with respect to a, we only need to 
prove the result for a = 1 (which we assume from now on). The recursive structure 
of the CMP is particularly suited for using renormalization arguments so it is not 
surprising that most of our proofs use a multiscale analysis. In the present case, we 
consider events of fhe form 

Z(n, 7 ) = ([[ 0 ; n [[) contains a cluster with a least jn elements} . 

Note that we consider here the CMP on the sub-graph [[0; n [[ which does not necessarily 
coincide with the restriction of “^(IN) to the interval [[0; n[[. We call these events good 
because of fhe following inequalify which holds whenever | < 7 < 1 : 

P{£( 2 n, 7 )}>l-{ 1 -P{£(n, 7 )})". 
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Figure 6: Red intervals represent bad events. Top: The bad event is on the interval [[?7„;2t?„[[ and 
Q cannot reach over to merge with Q. Bottom: The bad interval is inside [[2ti„;— 2t;„[[ and 
Q and Q merge together. 


To see that, apply the cluster merging procedure on [[0; n [[ and on [[n; In [[ separately. If 
each of these intervals contains a big cluster of size greater than ^n, then the rightmost 
vertex on the big cluster of the first interval is a distance smaller that | < 7 n of the 
leftmost vertex of the big cluster of the second interval. Therefore, these two clusters 
will merge together when performing the CMP on [[0; 2n[[ and give birth to a cluster of 
size at least 2'yn. 

Define = 20 • 2” for n > 0. We also define the sequence (7n)w>o by 7 o = 9/10 
and 7 „+i = 7 „ (^1 - . We can write 

n —1 ^ fi—1 2 ^ 2^1 9 1 8 

“ 10 ^ ~ - 10 ~ ~ W ~ - W ~ 10 ~ W' 

so that 3/4 < 7 „ < 1 for every n. Now, fix e > 0, by continuity, it is possible to chose 
p G (0,1) close enough to 1 such that P {£ (20,9/10)} > 1 — e. In the following lines, 
we will prove by induction that, for n > 0, 

P{£{Vn,7n)} > 1 - e/4” 


1 aer 

where Vn = uq • ui... Un- 

Assume that the inequality holds for some n > 0. The interval [[0; [[ is divided 

into Un+i sub-intervals of size Vn- Suppose that each of those intervals contains a big 
cluster of size at least except maybe one of them which we call the bad interval. As 
we already mentioned two clusters of size at least jn'Vn belonging to two neighboring 
intervals will merge since 7 ^ > 3/4. This means that all the big clusters belonging 
to sub-intervals on the left hand side of the bad interval will merge together into a 
cluster denoted by C;. The same thing happens for the clusters on the right hand 
side of the bad interval creating a big cluster C^. Now, the only case where C; and 
Cy do not merge together is when the influence radius of one of them cannot "reach 
over" the bad interval. This happens only if the bad interval is located at 2vn [[ 
or [[(wn+i — 2)vn', {un+i — (see Figure|^for an illustration). In any case, at least 
Uyi^i — 2 sub-intervals containing a big cluster merge together. Thus, the CMP on 
[[ 0 ; [[ contains a cluster of size at least 

{Un+l 2)7Mt’M = (1 J 'In’^n+l ~ T«+1^m+1- 

V Wfi+l/ 

Recalling that this occurs whenever there is strictly less than two bad intervals and 
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using the fact that the weights {r{x) ) are i.i.d., we get, using the recurrence hypothesis, 

1 - P{£(z;„+i, 7 „+i)} < P{Binomial(u„+i,l - P{£(z;„, 7 „)}) > 2} 

< uJ+i(l-P{f(u„,7n)})^ 

40^e^ £ 

< - < - 

— 4 W 4 n+l 

provided that e is chosen small enough. Thus, we have proved that, for e > 0 arbitrarily 
small, we can always find 0 < p < 1 such fhaf 

P{£(r'„, 7 „)} > 1 — £/4" foralln>0. 

We now infroduce fhe "anchored" good evenfs 

So{n, 7 ) = {in ^([[ 0 ; n[[) fhe clusfer confaining verfex 0 has af leasf jn elemenfs} . 

Using similar arguments as before, we see that the event £o(^w+i/7w+i) is realized 
whenever the following conditions are mef: 

• Soivn^Jn) is realized; 

• fhe inferval \vn',2vn [[ confams a clusfer of size af leasf ynVn', 

• For every 2 < k < Un+i except maybe one, the interval p:Vn', {k + 1)vn[[ contains a 
cluster of size af leasf ynVn- 

This implies 

P{£o(»ii+l.7n+l)} > P{f'o(!’ii. 7 n)} (1 - {f) P{Bmomial («„+!- 2 ,^) < 2 | 

> P{i-o(!>,„7„)}(l-{7) 

therefore, since 7 „ > 3/4 for all n, we gef fhaf 

liminf P {£^o(^w/3/4)} > 0 

n^oo 

which proves that, with positive probability, the cluster is infinite whenever p is 
sufficiently close to 1 . □ 

As we will see in the next section, for any a > 1, we have Pc{0i) > 0 for Bernoulli 
CMP on any d-dimensional lattice, d > 1. For the time being, we provide an alternative 
proof which works only in the one-dimensional case and for a = 1 but has the advantage 
of providing an explicif lower bound for fhe crifical percolafion paramefer. 

Proposition 3.5. Consider Bernoulli CMP on Z with a = 1. We have 

Pc(l) > 

The proof is based on fhe following dualify lemma. 
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Figure 7: Illustration of the proof of Proposition 2.12 In this configuration, 0 and L are both in a 


cluster composed of vertices of [[0;L]]. Vertices R and n +1 are also both in a cluster composed of 
vertices of |R;n + 1]]. In the next step of the cluster merging procedure, these two clusters merge. 


Lemma 3.6. We assume here that a = 1. Set In = |IO;n]] and fix a sequence of weights 
(r(x), X G In) ^ {0,1}^". Suppose that the CMP ^ {In, r) inside In is such that site 0 and site 
n belong to the same cluster. Then, the interval !„ (seen as a subset ofZ) is stable for the CMP 

constructed with the reversed weights f{x) = 1 — r{x). 


Proof of Lemma 3.6 The proof goes by induction. The case n = 1 is trivial. Suppose now 
that the result holds for every k < n. We consider In+i with a sequence of weights 
satisfying the assumptions of the lemma. Since sites 0 and n + 1 are in the same cluster, 
we have in particular r(0) = r{n + 1) = 1. According to the cluster merging procedure, 
the partition is obtained by successive application of the merging operator. We 

remark that there are at most n true merging needed to construct (because the 

number of clusters decreases by one after each true merging). Let (xj-, i/jt) be a sequence 
of pairs of sites describing a merging history for ^(7„_|_i) i.e. 


^{Ifl + l) — ^X,n,ym ° O • • • O Mx;^,yi(7„ + l) 

(where we identify In+i with its trivial partition). Let k < m denote the index where 
the operator Mxj^,yy. merges together the cluster containing 0 with the cluster containing 
n + 1 and let C be the partition of In+\ obtained just before this merging occurs. We set 


L A max Cq and R = min Cn+i 


i.e. L is the rightmost site of the cluster of C containing 0 and R is the leftmost site of 
the cluster containing n + 1. See Figure]^ for an illustration of this configuration. Since 
these two clusters can merge and since the weight of a cluster is equal to its number of 
sites (because each weight is either 0 or 1 ), we have 

R — L < min{L + l,n + 2 — R). (16) 


Furthermore, we observe that all the merging used to construct the cluster Cq containing 
0 and L are also valid when considering the CMP inside [[0; L]]. This means that the CMP 
on [[0; L]] has a cluster containing both 0 and L. Similarly, the CMP inside [[R; n + 1]] has 
a cluster containing both R and n + 1. Applying the induction hypothesis, we deduce 
that the subsets [[0;L]] and [[7?;n + 1]] are stable for the reversed weight sequence f. 
Moreover, the sum of the weights f of all the sites inside the center interval ]]L; 7?[[ is 
bounded above hy R — L — 1. Combining this with (16), we find that 


f{x) < min(L,n + 1 - R) < d(]]L;R[[ , Z \ [[0;n + 1]]) 

xe]L-,Rl 


and we conclude using Proposition 2.15 that [[0;n + 1]] is indeed stable set for the 
reversed weight sequence f. □ 
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Proof of Proposition |3.5| By symmetry, the probability of having an infinite cluster im- 
bounded towards +oo is the same as that of having an infinite cluster unbounded 
toward — oo. By translation invariance, these probability are either 0 or 1. Thus, if the 
infinite cluster exists, it is a.s. unbounded towards both +oo and —oo. 

We fix p = ^ and we suppose by contradiction that Pc ^ P- Thus, ^ contains an 
infinite cluster a.s. which we denote by C“. Let {xn,yn)n>o be a (random) sequence 
of pairs of vertices associated with the cluster merging procedure i.e. such that lim„ 

C” = ^ where 


(we identify Z with its finest partition). We can find an increasing function cp such that 
and i/(p(„) both belong to C“ for every n. The sequence of sets C” = increases 

weakly to as n tend to infinity. Let = min Cn and Rn = max C„ be the leftmost 
and rightmost vertices of C”. The cluster C“ being unbounded in both directions, 
and Rn diverge respectively go to — oo and +oo as n —> oo. Moreover, by construction L„ 
and Rn are in the same cluster if we consider the CMP inside [[L„; P„]]. Thus, according 
to Lemma 3.6 the intervals are stable sets for the reversed weight sequence. 

Moreover, the union of these sets exhausts Z so Theorem |2.20| states that the CMP on 
Z with weight sequence r contains only finite clusters. But since p = j, the sequences 
(r(x), X G Z) and (r (x), x G Z) have the same law which leads to a contradiction. □ 


Remark. We do not believe the lower bound 1/2 of Proposition |3.5| to give the critical value for 
Bernoulli CMP on Z. In fact, numerical simulations suggest that 


pc(l) ~ 0.65. 


We ask the question: does there exist an explicit formula for this critical parameter? 


3.2 Phase transition on for continuum CMP. 

In the previous section we have shown that, on any infinite graph, there exists a super¬ 
critical phase where an infinite cluster is present. In this section we prove that, on finite 
dimensional lattices, there also exists a sub-critical phase where every cluster is finite. 

Proposition 3.7. Let d > 1. Consider continuum CMP onG = "Zf with expansion exponent 
a > 1 and i.i.d. weights distributed as AZ. Suppose that E[Z^o] < +oofor /3o = (4xd)^. Then, 
we have 

Ac(x) > 0. 

By coupling, this proposition implies that the critical parameter for Bernoulli perco¬ 
lation is strictly positive in any dimension: 

Corollary 3.8. Consider Bernoulli CMP on with expansion exponent oc > 1. We have 

pM) > 0 . 

Remark. The moment condition is obviously not optimal. We conjecture that a sub-critical 
phase exists whenever Z admits a moment of order ad -\- efor e small enough. Conversely, ifZ 
does not have moments of order ad — e, then the maximum influence radius of sites at distance 
N from the origin is much larger than N and one can show that the CMP always contains an 
infinite cluster. 
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Figure 8: 


Left: An event S{N), the grey set is stable. Right: Constructing S{LN) from S{N). 


Proof of Proposition |3.7| The proof is again based on multiscale analysis. We define the 
"good event": 

S{N) = /there exists a stable set S such that [[N/5 ; 4N/5]]‘^ C S C [[1 ; N]]^ 


We will show that, when A is small enough, events £{N) occur for infinitely many N 
almost surely 

First, we remark that, since a union of stable sets is itself stable, for N, L C N, the 
event £{LN) can be constructed by translating boxes of side length N inside which 
£ (N) occurs. In order to cover the gaps on the boundaries of the smaller boxes, we 
cover the larger box using translations by N/2 for each coordinate, accounting for a 
total of (2L — lY boxes of volume N'^. See Figure illustration. Of course, not 

every one of these (2L — 1)'^ boxes covering [[1; LN \ needs to satisfy £ (N) in order for 
£ (LN) to be realized. We will say that the boxes for which £ (N) does not hold are "bad 
boxes". For k > 0, we define the event: 


4(N,L) 


at most k of the (2L — 1)'^ boxes of side length N 
covering [[1; LN^ are bad boxes 

We focus on events of the form £(Rn) and £k{Rn, T„+i) where 

Ln = 2"” and = Li... L„ 

with c = 2a.d + 1. We remark that 




We also define 


R — 2'-+'^ +...+c” — 2 "^ c-i—^ < J 


« o-2dc"+i 


(17) 


= 2 


The next lemma estimates the number of bad boxes among the boxes of side length Rn 
used to cover the larger box of side length Rn+i = Ln+iRn- 

Lemma3.9. Seiko = |'2^+^(c + 1)] . Suppose that for some n > l,we have 


F{£{Rn)}>l-en. 
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Then, it holds that 


P ^n+l) } — 1 2^”+!' 

Proof. The box is covered by (2L„+i — 1)^ smaller boxes of side length R„ 

shifted by multiples of Rn/2. These boxes intersect each other so the individual events 
that they are "good boxes" are not independent. However, we can partition this set of 
boxes into 2“^ groups such that each group contains at most mutually disjoint boxes 
{i.e. we put in the same group boxes translated by multiples of R„ in each direction). 
Therefore, if the event fails, one of these 2^ groups has to contain at least 

k' = \> 2c + 2 bad boxes. Using the fact that the events of being bad boxes are 
independent within each group of boxes, we find that 


l-P{4„(R»,C,+l)} < 2'*p|Binom(Lf+„£„) >*:'} 

^ ^2d+'^+dPc"+^-2dk'c'‘+^ 

2 

_ lnd+l-fc'dc"+U2dc”+2 ^ 

— 2^ ^”+1 

1 

^ 2^n+l 


where we used that d + 1 — k'dc^^^ + 2rfc”+^ < 0 when k' > 2c P 2 for the last inequality. 

□ 


We introduce another family of events meant to control to total sum of the influence 
radii inside a box: 

A{N,L)^^[ E r{x)<Ll] 
xell}Nf 

and we set 


g{N,L)^^ 


each of the (2L — 1)^ sub-boxes of side length N 
covering the box [[1; NL]]^ satisfies A{N, L) 


Lemma 3.10. There exists uq such that for all n > uq, uniformly on A < 1, we have 


P{g{Rn,Ln+^)}>l 


1 

2^n+l 


Proof The proof uses only crude estimates and imion bound. First, when A{N,L) fails, 
then there is at least one site x G [[1; such that r{x) > Consequently, we have 


1 


P {A{Rn/ Tfx+l) } 


< 


< 


< 


Rip{\z>T^ 

(Wl)^ 


E[Z^^]{Ln+i) 


c—1 
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where we used Markov's inequality and A < 1 for the second inequality and ( [T7l > for 
the last one. Now, using again the union bound, we find that 


l-P{g{Rn,L„+i)} < 

< 


(2L„+i-l)^(l-PM(R„,L„+i)}) 

2^E[Z/^o]e„+i(L„+i)‘^+^-^+2^A 


Recalling that c = 2a.d + 1, the exponent of Ln+i in the formula above is equal to 

3d+ — +4:0id^-^ < 8ad^-^ < 0 
2a. 2a 2a 

which completes the proof of the lemma. 


□ 


Lemma 3.11. Let ko be as in Lemma 3.9 There exists n\ > uq such that for all n > ni, 
uniformly in A < 1, it holds that 

if^nr Ln+l) n g{Rfi/ Ln+l) C E(^Rn+l)■ 

Proof The idea behind this inclusion is the following: on £i^^(Rn, L„^i), there are at 
most Lq sub-boxes which do not contain a stable set. On the other hand, on g{Rn, L^+i), 
the influence radius of these non stable sets is negligible compared to the diameter of 
the big box. Therefore, either a bad box is close to the boundary and it does not interfere 
with a stable set in the center of the box or it is at a macroscopic distance from the 


boundary in which case it is contained in a stable set according to Proposition 2.15 
More precisely, suppose that we are on the event £]^^{Rn, L„+i) (A g{Rn, L„+i). Let 
k < Rq he the number of bad boxes. We denote them by Bi ,..., and set 


R = 


E K^)- 

xeUiB, 


Let n > nQ. Since we are on g{Rn, Ln+i)> Jensen's inequality gives 

R“ < E ( E 

i=l \xeBi 

We consider a neighbourhood around each of the boxes of the form: 

Bi = {x G d{x,Bi) < R“}. 

From the triangle inequality, it comes that diam(B;) < 2R“ + Rn < 3 /CqL„_|_i. Conse¬ 
quently, for n large enough, 

f;diam(B;) < 3(:;+'L„+i < %i. (18) 

i=l 

We say that a bad box B, is cormected to the outside if there exists a path from some 
vertex of Bj to some vertex outside of the box [[R„+i/10; 9R„_|_i/10]]“ which stays 
inside \Jj Bj. 


j < KLn+l. 
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10 5 

Figure 9: In grey: the union of all the stable sets inside good sub-boxes. In red: the bad boxes B; 
and their neighborhoods B, which are connected to the outside. In blue: the bad boxes Bj and their 
neighbourhoods Bj which are not connected to the outside.The union of the grey set and the blue 
bad boxes is a stable set which contains the inner region delimited by the black square. 


Let So be the union of all the stable sets of the good sub-boxes of [[1; R„_|_i]]^. Let Si 
be the union of Sq together with all the bad boxes that are not connected to the outside. 
See Figure|^for an illustration. The set Sq is stable as a union of stable sets. Moreover, 
by definition W = Si \ Sq is included in the union of all the bad boxes hence 

E < R- (19) 

xew 


Since we added in Si only the bad boxes which are not connected to the outside, we 
have 


B(W,R“) C Si. 


( 20 ) 


In view of Proposition 2.15[ inequality ( [T^ combined with ( [20] ) implies that Si is also 
stable. Finally, inequality ^18|) says that any bad box connected to the outside is at 
distance at most from the outside of the box [[R„+i/10 ; 9P„_|_i/10]]'^. This implies 


SiD[[P„+i/5;4P„+i/5f 


which completes the proof of the lemma. □ 

We can now finish the proof of the proposition. As A tend to zero, we have Ar (x) —f 0 
for all X G Since 8 (Rni) depends only the values of r for finitely many sites, it 
follows that 

lim P{£ (R„.)} = 1 

A^O 

SO we can fix A > 0 such that 


R{8 (Rn,)} >l-en,. 

We prove by induction that the same inequality holds for all n > ni. Indeed, if it holds 
for n, then, combining Lemmas 3.9|3.10 and 3.11[ we find that 

F{8{Rn+l)} > l-{l-P{8t^{Rn,Ln+l)})-{l--P{G{Rn,Ln+l)}) > 1 - Sn+l- 
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By translation invariance and the Borel-Cantelli lemma we conclude that 
P I there exists a finite stable set containing [[—N;N]]‘^| = 1 
for every N which implies that the CMP has no infinite cluster. □ 

3.3 Phase transition on random geometric graphs and Delaunay tri¬ 
angulations 

We now explain how to extend Proposition |3^ for the degree-weighted CMP (Model 
1^ when the graph G is either a random geometric graph or a Delaunay triangulation 
constructed from a Poisson point process V in with Lebesgue intensity. First, we 
quickly recall the definition of these graphs: 

• Geometric graph with parameter R > 0. The vertex set is composed of the atoms of 
the point process V and, for any pair of points x,y ^ V, there is an edge between x and 
y if and only if ||x — i/|| < R, where || • || denotes the Euclidian norm in R'^. If R is above 
the critical parameter for continuum percolation, then this graph has a unique infinite 
connected component (see for instance mi). We assume this is the case and denote this 
graph Q{R,V). 

• Delaunay Triangulation. For any x ^ V, we define the Voronoi cell of x as the set of 
points of R'^ which are closer to x than to any other point of the Poisson point process: 

Vorp(x) ^^jzGR'^: ||x — z|| < ||i/ — z||Vi/g7^|. 

The Delaunay triangulation of V is the dual of the Voronoi tessellation: its vertex 
set is again the set of atoms of V and two vertices share an edge if and only if their 
corresponding cells are adjacent (i.e. they share a d — 1 dimensional face). We denote 
this graph by D(7^). Since the points of the Poisson process are almost surely in general 
position in R'^, this triangulation is also characterized by the following property: for 
any simplex of D('P), its circumscribed sphere contains no point of V in its interior. 

Proposition 3.12. Consider the CMP on Q(R, V) or on ViV) with weights given by 

r(x) deg(x)l^deg(^)>^j 

and with expansion exponent x > 1 Then, we have 

Ac(x) < 00. 

The proof is very similar to the proof of Proposition |3.7[ It would be redimdant to 
write everything again so we will simply point out the modifications needed to adapt 
the proof for these random graphs. The main difference (and difficulty) in our new 
setting comes from the fact that two portions of the graph included in disjoint domains 
of R'^ are not independent anymore. The modification needed for the geometric graph 
are minor since sub-graphs included in domains separated by a distance larger than R 
are still independent. The situation is a little more complex for Delaunay triangulations 
so we will concentrate on this case. Details for the geometric graph are left to the reader. 


28 




Sketch of the proof. Recall the parameters of the proof of Proposition|3.7[ 


Rn = Ll . . . Ln) 


r ^ « 9-2dc”+i 

— Z , tfi — Z 


for c suitably chosen and depending only on a and d. We adapt the definition of good 
events to the new setting: 


SiRn) = 

We want to prove by induction that for n large enough 


there exists a stable set S such that 
'V{V) n [R„/5,4R„/5]'^) CSC (v{V) n [0,Rn 



Once we have proved that the three lemmas 3.9 3.10 and 3.11 hold for our random 
graph, the result follows mutatis mutandis. 


Lemma 3.11| is a statement concerning stable sets of the CMP and it does not really 
depend on the particular nature of the graph. It is straightforward to adapt it to the case 
of Delaunay triangulafions and other random graphs embedded in R'^. 


Lemma 3.10 is also easily translated to our case. This lemma gives an estimate for the 


probability that every box has a reasonable total radius. As we already pointed out, it 
uses only the union bound so it does not require any kind of independence assumption. 
It also requires that the weights admit polynomial moments of high enough order. 
This is not a problem here since the typical distribution of the degree of a site for the 
geometric graph and the Delaunay triangulation has, in fact, exponential moments: this 
is straightforward for Q{R,V) since the degree of a site is bounded by the number of 
atoms of V inside a ball of radius R. For V{V), this result is proved in llTSi (and it also 


follows from Lemma 3.13 which we will prove later on). The last point to check is that 
we can upper bound the number of vertices inside a box. For the lattice case "Zf, this 
number was deterministic and equal to the volume of the box. In our new setting, it is 
still straightforward since the number of vertices is simply the number of atoms of V 
and hence if follows a Poisson distribution with parameter equal to the volume of the 
box. In particular, this distribution has light tails which is more than we need. 

Lemmarequires a bit more work. It controls the number of bad boxes in [0, R„+i]^. 
Recall that it states that, when 


P{£{Rn)}>l-en. 


then it holds that 


where 


£k{Rn> T 


n+1. 


R {£k{Rn/kn+l)} ^ 1 2^”+l 


at most k of the boxes of side length R„ 
covering [0,Rn+i]‘^ are bad boxes 


For the lattice Z‘^, we proved this by partitioning the set of small boxes covering 
[0, Rn+i]^ into groups containing disjoint boxes. Then we used the fact that good events 
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Figure 10: The blue circle circumscribes a triangle face {x,\j,z) T>{V)) hence it contains no 

point of V in it interior. The half disk H contains a red box of macroscopic size included in the 
annulus A. 


were independent for disjoint boxes to compare the number of bad boxes with a binomial 
distribution. We cannot do this directly now since we do not have independence of 
events occurring in disjoint boxes any more. This problem is easy to overcome for the 
random geometric graph as we simply partition the set of boxes into more groups in 
such way that two boxes in the same group are at distance at least R. The adaptation of 
the lemma in the case of the Delaimay triangulation is a bit more involved. Consider 
the events 




{ There are no edges linking vertices of 7^ n [f]N, (1 — ri)Ny 
with vertices of T’ n \ [0, N]^) 


Lemma 3.13. Let tj < 1/5. There is an event I{N,rj) depending only on the points ofV 
inside the annulus Atj = [0, N]^ \ [pN, (1 — p)NY with the following properties: 

i{N,p) C T{N,p), 

P|X(N,//)} > 

where C > 0 is a constant depending only on p and d. 

Let us first explain how we use this result to get the desired estimate. Let denote 
the probability that a small box does not satisfy I{Rn, 7 ) • By union bound, the probability 
that one of the boxes covering the larger box does not satisfies I{Rn,p) is smaller than 
d„(2L„_|_i — 1)'^ = o(e„_t_i). Thus, we condition on the event that every sub-boxes satisfy 
I{Rn,p). Then, the events when the small boxes are bad become independent for disjoint 
boxes and have a probability uniformly smaller than = e„ + ^ e„. Thus, we can 

again use a comparison with a binomial random variable with probability of success 
1 — ejj and the rest of the proof is the same as in the lattice case. □ 


Proof of Lemma 3.13 Take x e V H [pN, (1 — p)Nf and y ^ V D (R“ \ [0, N]“) and 


suppose that there is an edge between x and y in V{V). There exists an (hyper)sphere S 
circumscribing a simplex of the triangulation and having the segment [x, y\ has a chord. 
By definition of the Delaunay triangulation, this sphere has no points of V in its interior. 
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Let [^',y'] = [x,i/] n Arj. The sphere S' having the segment [^',y'] as a diameter has an 
hemisphere H included in S. It is now easy to convince oneself that there is always a 
square box B of side length SN included in the intersection A n H of the annulus A and 
the hemisphere H {c.f. figure]^. By construction this box does not contain any point of 
V. Moreover, S depends on y and the dimension d but not on N. Thus, we just need to 
construct an event I for the point process V such that any possible square box of side 
length SN inside A^ contains at least one atom of V. This is achieved by partitioning 

the annulus into boxes of side length yN/k where k = \2y/S^ and then requesting that 
the point process has at least one atom in the interior of every box. Again, the number 
K of boxes needed depend only on y and d. The probability that the event I constructed 
in this way fails is the probability that there exists an empty box. By union bound, this 
is smaller than K exp { — (yN/kY) and the lemma follows. □ 

Remark. 

• The multiscale technique used to prove the existence of stable sets and therefore a sub- 
critical phase is quite robust. Here, we only used it for three particular weighted graphs 
but it is easy to convince oneself that it can be applied for many other graphs that can be 
embedded into a finite dimensional space in a "nice way" (for graphs with exponential 
growth, using a multiscale method seems much more challenging). For example, the 
previous arguments work for general point processes provided that the intensity is bounded 
away from 0 and oo. We can also consider random radii when constructing the random 
geometric graph as long as the distribution has very light tails. 

• Another family of graphs that we think would be interesting to study are the infinite 
uniform planar maps (such as the uniform infinite planar triangulation). These are graphs 
for which the degree of a typical site has exponential tails so we expect again that > 0 
for degree weighted CMP. However, this will require more work and it is not clear (to us) 
what embedding into should be chosen in order to use a multiscale argument. 


4 Connection with the contact process 

In this section, we make rigorous the heuristic given in the introduction by relating the 
almost sure extinction of the contact process on an infinite graph to the existence of a 
sub-critical phase for degree-weighted CMP on the same graph. 


4.1 The contact process on a locally finite graph 

Recall that G = (P, E) is a locally finite and connected graph. Fix a parameter A > 0 
which we call "infection rate". The contact process ^ = (^(t), t > 0) on G is a continuous 
time Markov process taking values in {0,1}^ with transition rates given, for each x C V 
and A C P, by 

A ^ A — {x} at rate 1, (21a) 

A A U {x} at rate Ally G A, d(x,i/) = 1}|. (21b) 


Thus ^(f) represents the subset of infected sites at time t and ( |21a| > means that each 
infected site recovers at rate 1 whereas ( |21b[ ) states that each site, while infected, emits 
independent infection vectors along its adjacent edges at rate A. 
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When the graph G has bounded degree (in particular when it is finite), classical 
theorems concerning interacting particle systems show that these transition rates define 
a unique Feller process (see for instance [Zl or |j8l for details). However, in our setting, 
the graph G usually has unbounded degrees and we need to be a bit more careful 
when defining the contact process. In order to do so, we use the classical "graphical 
construction" which we briefly recall, see for example [jZl, p32 for additional details 
about this representation. 

For each x ^ V, let Nx denote a Poisson point process with intensity 1 on R-|-. For 
each oriented edge {x,y), let Nx,y denote a Poisson point process on ]R+ with intensity 
A. We assume that all these Poisson processes are independent. Consider H = V x R+. 
For each x G R, put "recovery" marks on the time-lines {x} x R_|_ at the position of the 
atoms of Nx- For each oriented edge (x,i/), put arrows from (x, t) to (y, t) at the times t 
corresponding to atoms of Nx,y. Following Liggett BZl, we call active path a cormected 
oriented path in H which moves along the time lines in the increasing t direction, jumps 
from a site to a neighbouring one using the oriented arrows but never crosses any 
recovery mark. Then, we define the contact process (^(f), f > 0) on G starting from an 
initial infected configuration ^(0) = A in the following way: 

^(t) = {x G V, there exists a finite active path from (y,0) to (x, t) for some y G A} 

( 22 ) 

(see Figure]^ for an illustration of this construction). This construction defines the 
contact process for all time t > 0 in terms of a particular oriented percolation process on 
G X R+. Let us point out that without any additional assumption on G, it is possible that 
the process starting from a finite number of infected sites blows-up {i.e. creates infinitely 
many particles) in finite time (this corresponds to having an infinite percolation cluster 
in a slice G x [0, f]). However, we will not be concerned with this case as it will be ruled 
out by the additional assumptions that we shall make on the graph G. 

We use the notation to emphasize the starting configuration A of the process. We 
will also need to consider the contact process defined on a subset W C R of the vertices 
which we will denote by ^ This process is constructed using the same graphical 
representation and by keeping only the infection arrows linking vertices inside W. 
Finally, we define the number of infections exiting W up to time T as the number of 
oriented arrows in the graphical representation of the form (x, f) —> (y, f) with x G W, 
y ^ W and t < T that are reached by an active path starting from a vertex in ^(0) and 
which stays inside W. 

A nice property of the graphical construction is that it provides a natural coupling 
between processes defined in distinct sub-graphs and/or with distinct initial infected 
sets. More precisely, for any W C and any A C A', we have 

^|Jv(0 C for all t > 0. (23) 

If (W„) is an increasing sequence of finite sub-graphs of V such that lim j" W„ = V, it 
follows from this coupling and ( [2^ that 

^(f) = lim t ^|w„(0 for all t > 0. 

This means that the process ^ defined on the infinite graph corresponds to the weak 
limit of the contact process defined on any increasing sequence of finite sub-graphs W„. 
In fact, the main theorem of this section is stated in terms of ^ on the infinite graph 
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Figure 11: Graphical construction of the contact process on Z. Green crosses are recovery times, 
red arrows are infection times. The blue graph represents the history of the contact process started 
with site 0 infected: we can see that at the time t only —1, 1, 3 and 4 are infected and that it dies 
out in finite time. 


but we could restate it in terms of the contact processes ^ restricted to finite subsets 
which are, ultimately, the only processes we will consider during the proof. 

The graphical construction of the contact process gives a direct proof of the important 
self-duality property of the model: since the Poisson processes (Nx) and {Nx,y) are 
invariant by time reversal, it follows that, for any fixed time f > 0 and any two sets A, B, 
we have 

n B ^ 0} = P{^S(f) n A 0}. (24) 

We can now state the main theorem of this section which provides a sufficient 
condition on the geometry of a graph G to ensure the existence of a sub-critical phase 
for the contact process. 

Theorem 4.1. Let G = {V,E) be a locally finite connected graph. Consider the degree-weighted 
CMP on G i.e. with weights given by 

r{x) =%(^)l{deg(x)>A}- (25) 

Suppose that for some expansion exponent a > ^ nnd some A > 0, the partition ^{r,ci) has 
no infinite cluster. Then, the contact process on G has a sub-critical phase: there exists Ag > 0 
such that, for any infection parameter A < Ag, the process starting from a finite configuration of 
infected sites dies out almost surely. 

Remark. Let us make a few comments on Theorem ^l] 

• First, we find it remarkable that, in a way, all the geometry of the graph needed to prove 
the existence of a sub-critical phase is encoded in the merging procedure: the radii r{x) 
give the degrees sites but provide no information on the local shape or growth of the graph 
around a site. In particular, the theorem requires no assumption on the growth rate ofG. 
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• The exponent 5/2 is not optimal. However, the proof we describe cannot yield an exponent 
smaller than 2 so we did not find it worth the effort to clutter the proof with additional 
technical details for very little gain. In order to get an exponent close to 1, we believe that 
one needs a better understanding of the inner structure of clusters. The real challenge is to 
prove (or disprove) the theorem for oc = 1. 


Theorem 1.1| stated in the introduction of the paper is now a consequence of the 


Theorem above and Proposition 4.1 


Corollary 4.2. The contact process on a random geometric graph or on a Delaunay triangulation 
admits a sub-critical phase. 


Let us give a rough description of the proof of Theorem |4.1[ The basic idea behind 
the theorem is that when the infection parameter A is very small, the sets where the 
contact process is locally super-critical are the big clusters of the CMP. Yet, when the 
CMP has no infinite cluster, these big clusters look like islands surrounded by an ocean 
of small degree sites; and, on this ocean, the contact process is sub-critical and dies 
out quickly. We will prove that we can find a neighbourhood S aroimd each cluster C 
that will compensate the super-critical activity inside the cluster. More precisely, we 
will show that when an infection reaches the big cluster C, even though it will likely 
generate many infections before the whole cluster recovers, only very few infections 
will exit the neighbourhood S (less than one in average). Then, we can couple our 
process with a sub-critical branching Markov chain to conclude that the process dies 
out. The difficulty here is that we need these estimates to hold for every single cluster. 
Otherwise the coupling is useless since the branching process can survive locally on the 
finite graphs where the estimates fail. 

A natural candidate for the neighbourhood S is the stabiliser of the cluster. It turns 
out that we need to consider a slight modification of these sets in order to have more 
control on their size compared to the size of the cluster C. This is the purpose of 
Section where we define the notion of -stabilisers. We also prove in this section 
the key Proposition |4.5| which tells us that, indeed, these -stabilisers are large enough 
to dissipate most of the infections generated by their cluster. This is where we require 
Oi > 2.5 in order to have enough room to bootstrap the result from smaller //-stabilisers 
to larger ones. Again, the proof makes heavy use of the multi-scale structure of the 
CMP since stabilisers are themselves composed of smaller stabilisers. 

In Section 4.3 we introduce the particular branching process that we will couple 
with the contact process and present estimates for the extinction time and number of 
particles created that will be needed for the last steps of the proof. 

Finally, in Section 4.4 we put everything together, prove the main estimates and 
complete the proof of the theorem. 


4.2 //-stabilisers and the graph 


Recall that, according to Corollary 2.19[ the stabiliser Su of a subset U is equal to the 
union of all the clusters intersecting U together with all their descendants in the oriented 
graph on the set of clusters Fix 0 < tj < 1 and consider another adjacency relation 
on the set of clusters ^ given by 


C 4 C' 


C ^C' and d(C,C') < y(r(C))‘^. 
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V = 1/2 


i ^ tj = 1/S 

Figure 12: Examples of ?/-stabilisers. The underlying graph is the same as in Figure]^ Initial weights 
are in grey and weights of clusters are colored. Some ?/-stabilisers of the red cluster of total weight 8 
are displayed in red brackets. Notice that the 1/8-stabiliser is not a connected graph. 


For f] = 1, this corresponds to the previous definition and for t] < 1, it is a more 

restrictive condition so the oriented graph (^, • —/ •) is a sub-graph of (^, • —)• •) 
defined in Section [2^ In particular, it does not contain any cycle or any infinite oriented 
path. Mimicking the definition of stabilisers, we introduce: 

Definition 4.3. For any subset W C V, we call rj-stabiliser ofW the union of all the clusters 

fj 

of^ that intersect W together with all their descendants in the oriented graph (^, 
denote this set by S^, and write for 

Remark. Contrarily to 1-stabilisers, rj-stabilisers are not necessarily connected sets. See Figure 
[T2j/br an example. However, by construction, we still have the property that any two rj-stabilisers 
are either disjoint or one of them is included in the other one. 

We use //-stabilisers to define yet another new oriented graph, = {V,- —/ •), with 
the same vertex set V as the original graph G but with adjacency relation • —/ • given by 

x^y d(Sx,y) = 1 (26) 

i.e. the descendants of a site x in this new oriented graph are exactly the sites on the 
outer boundary of its //-stabiliser. Notice that if ^ has no infinite cluster, then every 
stabiliser is finite so the out-degrees in are finite. 

Remark. The relation —/ is defined on the vertex set V whereas —/ and —/ are defined on the 
set of clusters 

One of the main advantages of dealing with //-stabilisers instead of 1-stabilisers is 
that we have a precise control of their size which, in turns, provides sharp estimates for 
the distance between two adjacent sites in G'T 

Proposition 4.4. Let x,y such that x ^ y, we have 

yrf^xT < d{x,y) < 1 + ^rf^xT 


where 



—+J- 

2 “ - 2 1 - 7 



Remark. For // = 1, the proposition fails and the upper bound is of order rf/F’x Y'^^. 


(27) 
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Proof. By definition, any site z G V with d{^x,z) < f]r{^xY belongs to the //-stabiliser 
Sx- This proves the lower bound. For the upper bound, fix z ^ S^. There exists a chain 

of clusters = Cq —^ Ci —> C 2 —^ ... —^ C„ = for some n > 0. By definition, we 
have d(Q, Q+i) < rjr{CiY. This implies that r(Q_|_i)^ < rjr{CiY since otherwise Q and 
C;_|_i would have merged together during the CMP. Thus, for all 0 < i < n, 

r(CiY < vH^xY- 

In view of Proposition |2.6[ we find that 

diam(Q) < j^^^zr^r{^xY 

and the triangle equality yields 

n 

d{x,z) < diam(^x) ^ (d(Cf_i,C/) diam(C;)) 

i=l 

< + E YriV.r + 

;=1 

= 0 + 2 =^)) 
which yields the upper bound. □ 

Proposition 4.5. Assume that a = 2.5 and tj = 0.1. Define 

f{C) = r(C) 2 for any cluster C (28) 

Fix C G ^ and consider a chain 

xq ^ xi ^ ^ Xn where Xi ^ S^ \ Cfor aZZ z = 0,1,... , n — 1 

(we do not impose any restriction on which may either be inside Sq\C or on its outer 
boundary). We have, for any f G [1,a:], 

E > d{xo,Xn)^. (29) 

i=0 

Moreover assuming that d(xQ, C) = 1, d(x„, S^) = 1 and r(C) > 100, we have the stronger 
inequality, valid for 1 < /3 < 1.01, 

E !■(*■.,(30) 

i=0 

Let us give an interpretation for this proposition. Imagine that, for each vertex x ^ V, 
we must pay a price ffrfx)^ in order to travel along one of its outgoing edges in . 
Inequality ( |^ tells us that, in order to travel from some site x to some other site y, we 
must pay a total price at least d(x,yY^^- The second inequality ( [^ says that when /3 is 
close to 1 and when the cluster C is big enough, if we want to exit starting from a 
some boundary point of C, then any road staying inside <S^\C will cost more than the 
price f (C)^ required to travel directly from a site of C to the outer boundary of S^. 

The proof of the proposition is based on the concavity formula stated in the next 
lemma. 
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Figure 13: A cluster C and its ?^-stabilizer with a path {xq, ■ ■ ■,x„) in O’? as in Proposition 

x- 


In this example, the endpoint x„ in on the outer boundary of 


4.5 


Lemma 4.6. Lete,6 G [0,1]. We have 

inf I (x® + ... + x^) : G N, Xi + ... + x^ = 1, 0 < Xf < e| = g(e,0) 

whereg{e,6) = + (1 - > 1- 


Proof of Lemma 4.6 Set n = [Ij. Choosing = n + 1, xi = ... = x„ = e and Xn+i = 
1 — ne we see that the infimum is indeed smaller that ne® + (1 — ne)^ = g{£, 6). For the 
converse inequality, we use the fact that, for a > b, we have (a + x)® + (& — x)^ > a^ + b^ 
for any x > 0 and work by induction to transfer mass e on the first n values Xi,..., x„ 
and the remaining mass 1 — ne on x„_|_i. This argument also shows that g(e, d) > 1 for 
any choice e, 6 G [0,1] . The details are left out to the reader. □ 

Proof of Proposition |4.5| Set 

n—l r(W 1 “ 

Pt = YL and e = rnax 

i=0 


0<i<n-l D 


On the one hand, using Lemma 4.6 we find that 


i; = E 

i=0 i=0 




On the other hand. Proposition |4.4| states that rf(X;, x^+i) < 1 + 'pTf^XiY ^ 
hence 


1 


n—l 


D> -YL > -d{xQ,Xn] 

O' i=Q T 


( 31 ) 


Combining these inequalities, we find that 


E r(v.f > g 1 
;=0 
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4(XQ, Xyf^ ^ . 


1 

7 














Inequality now follows from the fact that g{e, ^/oc) >1 and 7 ~ 0.41 < 1. 

We now turn our attention to the second inequality. Under the assumption that 
d(xo,C) = 1 and d{xn,S^Q) = 1, the lower bound of Proposition 4.4 gives d{xQ,Xn) > 
rjr{CY — 1 hence 


n—1 


z =0 




= ^ e. 


V 1 

7 'yr{C)^ 


•(C)'*. 


Since Xq,...,X n-i are all in \ C, we have r{^XiY ^ t]t{CY and therefore 


but, according to 


ri'^x.Y < (2 + ?/^r(C)) 


D > ^d{xQ,Xn) > ^{r}r{CY — !)• 


These last two inequalities combined together yield 


e < 


2 + r]otr{C) 


— 1 
7 


Uj/r(C)“-l) 


= 7 


rj«-r{C) 


+ 1 


1 - 


rjr[C)‘' 


(32) 


For our particular choice of the parameters a = 2.5, r] = 0.1 and when r(C) > 100, one 
can check that e < 0.47 which, in turn, implies that for 1 < /3 < 1.01 


^ e. 


f V 


' X J \7 'yricy 


> 1.01 > y. 


This inequality together with (|32|) completes the proof of the proposition. 


□ 


4.3 A time-delayed branching Markov chain 

We now introduce another stochastic process that we will use to dominate the contact 
process. This process will be defined on the oriented graph but we give here a 
general description. Let H = (W, • —> •) denote an oriented graph. We assume that 
every site has finite out-degree. For each x G W, consider a real-valued random variable 

Tx > 0 

together with a family of integer-valued random variables 

Px,i/i '/■■•'/ ^x,yn ^ [[ 0 / [[ 

where { 1 / 1 ,... i/n} is the set of neighbor sites x —)• yi. We do not assume any indepen¬ 
dence between the random variables (tx, Bx,y ^,..., Bx,y„)- We call Time-Delayed Branch¬ 
ing Markov Chain (TDBMC) a continuous-time system of particles (Xt(x),x G W) such 
that: 
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• Xt{x) G {0,1,...} represents the number of particles at site x, at time t. We start 
from an initial configuration of particles (Xo(x), x G W) at time 0. Note that there 
may be more than one particle per site. 

• Each particle evolves independently of the others: when a particle is created at 
some time t, its survival time and progeny is independent of the evolution of all 
the other particles in the system at time t. 

• When a particle is created at some site x, it stays there for a random time fx after 
which it disappears while giving birth to Bx,y new particles at each neighboring 
site X — > y with 


fx. Bx,y„ Bx,yJ {Tx. Bx,y„ Bx, 


'f}/n ) * 


In order to define such a process at every time t >0, we must ensure that there can be 
no explosion in finite time. This is the case as soon as 

inf E[tx] > 0 (33) 

xew 


which will be a standing assumption from now on. Let us remark that if = 1 a.s. 
for all X, then X is a classical discrete time Branching Markov Chain. Another special 
case is when all the Tx's have exponential distribution; then the TDBMC is a continuous 
time Markov process. Notice however that despite its name, the process X does not in 
general satisfy the Markov property. 

We use the notation Ex[-] to denote the expectation for the process started at time 0 
from one single particle located at site x. Define also 


bx = E 


Lb,,, 

x^y 


, def E[Bx,y 

Ky = 


and 


def I r 

Ux = E T, 


X\r 


(34) 


with the convention \x,y = 0 if = 0. The next proposition collects properties of this 
process that we will use. 

Proposition 4.7. Let xq G W. Consider the TDBMC X started from a single particle located 
at site Xq. 

1. Let N be the total number of particles born in the TDBMC up to time +oo, we have 

00 / n—1 \ 

Exo[N] < sup 

w=o {xQ,...,Xn)e'P^Q V*=o / 

where 7^”^ is the set of all oriented paths of length n starting at xq. By convention, the 
product over an empty index set is equal to 1. 

2. Let T G [0, +oo] denote the extinction time ofX. We have 

00 / n—1 \ 

Exo[r] < E sup [uxnYl^^,]- 

n=0 {xo,...,x„)e'P^Q \ i=0 ) 
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3. Let U (Z W with xq i U. Consider a modification of the process where all the particles 
entering U are frozen (i.e. when a particle reaches U, it does not reproduce and it stays 
there forever). Then, starting from one particle located at xq, we have 


E 


*0 


Total number of particles frozen 
in U up to time t = +oo 


< sup (n&x,) 

(xo,...,Xn)ePTg \i=0 ) 


where is the set of all finite oriented paths {xq, ..., x„) starting at xq, with xi C W\U 
for i < n, and with G U. 

Proof Starting from one particle at site Xq and conditioning on its progeny, we get the 
relation 


E,.[N] = 1+ E E[B,„,.]E„[N] =1+ Y, i>.oA.v.,v,E„(]V]. 

Xq^Xi Xq^Xi 

Expanding this induction relation we get by monotone convergence 

00 n — 1 

n=Q {xo,...,Xn)eVSQ i =0 

00 / n —1 \ / n —1 \ 

^ E sup n J E n • 

n=0 y(xo,...,x„)€P?Q i=0 j y(xo,...,x„)6P?g !=0 j 

Recalling that {?^x,y) is a (possibly defective) transition kernel i.e. Yly^x ^x,y < 1/ the 
sum after the supremum on the right hand side of the last inequality is bounded above 
by 1 which completes the proof of the first statement. 

The proof of Item is obtained similarly starting from the inequality 

ExolT] < E[Txg] + YYi ^[^xo,Xi]'Exi[T] = Uxq + &xgAxo,xiExJr] 

Xq^Xi Xo^Xl 

which gives, using the same bounds, 

00 n — 1 00 / n —1 \ 

Exoi^"] ^ E E n E sup n ■ 

n=0 (xQ,...,x„)e'P^Q t=0 n=0 (xQ,...,Xn)e'P'f \ i=0 J 

We now prove ltem|^ The number of particles frozen does not depend on the Tx's, so 
we just need to consider the case Tx = 1 a.s. for all x C W. But in this case, the TDBMC 
is simply a discrete-time branching Markov chain. Let X(;(x) denote the number of 
particles at time n and site x for the discrete-time branching Markov chain obtained by 
freezing particles in U. Its transition kernel {Bx^y)x,yeW is given by 



l|x=i/} if G IT, 
Bx,y otherwise 


(here, we implicitly added a loop at each site x C U so that x —^ x). Now, define 


b 


/ 

X 


1 if X G Lf, 
bx otherwise. 


and 


„ /E[By/(,; 
*•!' \o 


if b'x > 0, 
otherwise. 
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Once again, we have a recurrence relation, namely 


y^x 


which implies 


n—1 


E»„ix;(x)]= E 

{xQ,...,Xn)e'P^^ i=0 


(the paths in above are considered for the new graph where there is a loop at each 
site of U). Summing over all x G Lf, we find that 




Total number of particles 
frozen in U up to time n 


+1 


< 


E n >’'xA„,. 

{xo,...,Xn)ePll^ !=0 

Xnell 

( \ ( 

k-1 

sup rr br- 

{xo,...,Xk)eV^, i=o ‘ 

\ k<n / 


E ff 

ixo,...,Xn)eP^^ t=0 

\ XnElI 


\ 




k-1 

< sup TT br-- 

{xo,...,Xk)eP^, i=o 

k<n 


We conclude the proof by letting n go to infinity 


□ 


4.4 Proof of Theorem I4.ll 

We can now state the main estimates which assert that, starting from a completely 
infected cluster C, the expected number of infections exiting the //-stabiliser for 
the contact process inside decreases faster than exponentially with respect to the 
cluster's weight. 

Proposition 4.8 (main estimates). Fix t] = 0.1 and x = 2.5. Suppose that the CMP with 
parameters given by ( [25l > has no infinite cluster for some A > 0. Then, there exists Ag > 0 
depending only on A such that, for any cluster C and any infection rate A < Ag, we have 


Total number of infections exiting 
for the contact process 


“ 2 


1.01 


(where f(C) = r(C) +2 as in ([28l>) and 


Extinction time of the contact process 


< e 


3r(C) 


(35) 


(36) 


The strength of these estimates is that Ag only depends on the geometry of G through 
a and A. Thus, the proposition gives bounds that are uniform for any cluster of any 
graph whose associated CMP has no infinite cluster. Let us first show how these 
estimates easily imply Theorem |4.1[ 
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Proof of Theorem |4. 1 1 By monotonicity of the CMP w.r.t. the expansion exponent a, we 

only need to prove the theorem for cc = 2.5. Recall that O'? = {V,■ —?■ •) denotes the 

oriented graph with vertex set V where the sites y such that x —> y are exactly those on 
the outer boundary of the y-stabiliser of x. Consider the following modification of the 
contact process where there may be more than one infection at each site. 

• We start from an initial infected vertex xq. 

• At time 0, we instantaneously infect every site of the cluster Then, we run a 
contact process inside Sxg, freezing every infection exiting the y-stabiliser. 

• At the time when the contact process inside Sxg dies out, for each infection that 
exited Sxg, we restart an independent process from the endpoint of the infection. 
This means that, for each infection with endpoint, say z, we instantaneously infect 
every site of arid then run an independent contact process inside Sz, freezing 
the all the infections exiting the y-stabiliser. 

• We construct the process for all times (or until extinction) by iterating this proce¬ 
dure. 


The freezing of particles in this modified process prevents us from coupling it at deter¬ 
ministic times with the real contact process started from the same initial infected site 
Xq. Yet, we can still construct both processes on the same probability space in such way 
that: 

1. The extinction time of the modified process is larger than or equal to the extinction 
time of the contact process. 

2. For any directed edge of the graph G, the total number of infections sent through 
the edge by the contact process is smaller than or equal to the number of infections 
sent by the modified process. 

This coupling can easily be achieved by using the same sequences of clocks on sites and 
oriented arrows for both processes. However, contrarily to the graphical construction, 
in this case, the time on a site (resp. oriented edge) runs only when the site (resp. start 
vertex) is infected. This ensures that the modified process will never miss any infection 
clock that the contact process uses. We leave the details to the reader. 


Looking at infections exiting y-stabilisers, we see that this modified process naturally 
defines a TDBMC X on the graph with transition kernel given by (with the notation 
of Section [43l >: 

Tx = Extinction time of the contact process (37) 

and for every y G R such that x —^ y. 


B 



total number of infections reaching 
site y for the contact process f,. L. 

\^x 


(38) 
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From the coupling with the contact process, we see that condition]^ on extinction times 
means that 


inf(t > 0, = 0) < inf(t > 0,Xt{x) = 0 for all x G V). 


Using the notation of ( |34| >, the main estimates translate to 

1 


bx < 

2 




(39) 


Therefore, in view of Item of Proposition |4.7[ we find that the expected total number 
of particles created in Z starting from xq is bounded by 


E 


n—1 


sup 


n=0 {xo,-,x„)eV^g z=0 


n < E • 


n=0 


(here, denotes the set of paths in starting from xq with length n). This means 
that X creates only finitely many particles a.s., hence its extinction time T is also finite. 
This in turn implies that the contact process dies out almost surely □ 


Proof of Proposition \4.8\ The proof works by induction on the weights r(C) of clusters. 
From now on, we fix 

a = 2.5, A > 1, f] = 0.1 

and we assume that the CMP defined by ( [25l > has no infinite cluster. Let us remark that, 
since, for every cluster C, the stabiliser is a finite set, the random variables inside the 


expectations in ( |35| ) and ( |36| > have exponential tails and are therefore finite. This follows 
easily from the graphical construction described previously. 

For every R > 1, there exist only finitely many graphs isomorphic to some t]- 
stabiliser where C is a cluster with r(C) < R. This see this, notice that r(C) < R 
implies that every site in the stabiliser has degree at most X V A. The number of sites in 
the cluster C is also bounded by R. Using Proposition |4.4[ we deduce that the number 
of sites in the stabiliser is also bounded by f{R) for some function / growing fast 
enough. This proves our assertion since there are only finitely many non-isomorphic 
graphs with at most f{R) vertices with degrees bounded by X V A. 

For each cluster C, when the infection parameter A of the contact process goes to 0, 
the expectations in ( [^ and ( [^ tend to 0 by dominated convergence. This means that 
we only need to prove Proposition 4.8 for clusters satisfying r(C) > X where X = X(A) 
may be chosen arbitrarily large. 

From now on, let X be large and A > 0 such that Proposition |4.8| holds for every 
cluster C with r(C) < X. We will prove that the same result holds, in fact, for every 
cluster C with r(C) < ^X and the result will follow by induction. 

Fix C G ^ such that r(C) < ^X. We use following the notation for conciseness: 


S = Si 
D = Sl\C. 


We decompose the proof in six steps. 
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1.01 


Step 1. For any xq ^ D with d{xo, C) = 1, we have 

Total number of infections exiting through dS 
for the contact process 

(remark that we do not count in this expectation the infections going from D into C) and 


< g-l.Olf(C) 


Extinction time of the contact process ^ 


XQ 

D 


< 2c 


(40) 


(41) 


where c is a universal constant. 


We consider again the TDBMC X on with transition kernel given by (371 and 

_ ]), where we freeze all the particles exiting D. Using the same arguments as in the 

proof of Theorem |4.1[ we can couple X started with one particle at Xq with the contact 


process in such a way that: 

1. The total number of particles in X frozen on the outer boundary of S is larger than 

^0 

D' 


the total number of infections sent outside of S by the contact process 


2. The extinction time of X is larger than the extinction time of the contact process 




XQ 

D- 


Now, since every cluster C' inside D is such that r{C') < rjr{C) < R, we can use the 

main estimate to upper bound the quantities {bx,x G D). Denoting by the set of 
paths in which start from xq and such that x, G D for i < n and G U \ S, we find 
with the help of Item of Proposition 4.7 that 

Total number of infections exiting through dS 


rXo 


for the contact process ^ 


< 


sup 

{xo,...,Xn)e'Pf^^ 



< 


sup 

(xo,...,Xn)e'Pf^^ 


2/1 ^ 


X,') 


1.01 


g-l.Olr(C) 


1.01 


where we used Proposition |4.5 for the last inequality This completes the proof of 
(|4^. The proof of the second inequality is similar. Let 7^”^ denote the set of paths in 
G't staying inside D, starting at xq and of length n. We use Item of Proposition |4.7| 
combined with Proposition |4^ to get that 

Extinction time of the! 

< sup 

n=0 ixo,...,Xn)e'P^Q 


Xq 


contact process 


-d{xo,Xn 


1.01 


< E 

n=0 {xo,...,x„)e'P^^ 

00 

< E 

n=0 (xo,...,x„)€T^g 
00 

< E 4 < 2c 

— nn — 

where c is the overall supremum on [0, oo) of the function x —7'3(l + x)a — x^. 



1.01 
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Step 2. We have 


P 


The contact process never 
sends any infection into C 


1 

- 2''(C)‘ 


(42) 


Let ei,... ,em denote the set of edges connecting D to C. Using the fact that the 
contact process has positive correlation at all times (c./. Theorem B17, p9 of [ 8 J) it is easy 
to check that 




sends any infection into C 



y vxix^w-^xi x^x xxxx p ^ /1 i- 

The contact process never 1 
sends any infection through e,. J 


Fix Ci = {xi —)• Zi) with X; G D and Z; G C. Let f denote the total time site Xi spends 
infected: 

poo 

Using the self-duality property of the contact process we find that 


E[/i] 


If F n {xi} ^ it 
If F [pl^it) n {D} ^ it 




dt 


E 


Extinction time of the contact process 


< 2c 


where we used Step 1 for the last inequality. Now, recall that an infection propagates 
through Ci whenever site Xf is infected and a clock attached to the oriented edge X; —> Zf 
rings. Since these clocks are independent of the contact process the expected 

number of infections exiting through Cf is bounded by 2cA. We can without loss of 
generality assume that 2cA < | in which case we have 

{ The contact process ^ 1 ^ never 1 ^ l 
sends any infection through Cj J 2 


The claim follows noticing that the number m of edges between C and D is bounded by 
r(C). 

Step 3. We have 


The contact process never 
sends any infection into D 


> 


2r{C)' 


and 


Extinction time off,' 


< 2''(C). 


(43) 
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This step is easy. We use a very crude estimate: the probability that no infection 
ever escapes C and that the contact process dies out before time 1 is larger than the 
probability that every (recovery) clock attached to a vertex of C rings before time 1 and 
no (infection) clock attached to an oriented edge with a start vertex in C rings before 
time 1. Therefore, for A < 0.1, using again that r(C) upper bounds the number of 
outgoing edges and vertices in C, we get 

f The contact process never sends any 1 \ '"(c) 1 

infection into D and dies out before time 1 ^ ^ 


Now, comparing the contact process with the modified process obtained by re¬ 
infecting every site of C at each integer time when there is at least one infected site, 
it follows from the previous inequality that the extinction time of is stochastically 
dominated by a geometric random variable with parameter This gives ( [4^ . 

Step 4. Denote by the total number of infections traveling through an oriented edge from 

C to D for the contact process started from the configuration A and restricted to S. Similarly, 
define ^ the number of infections traveling through an oriented edge from D to C and 

t-A ^ a a 
^ \S ~ ^ \S' ^ \s 

the total number of infections travelling between C and D. For any initial infected set A C S, 
we have ^ 

P{rfs>k} < forallk = 0,l,... (44) 

By stochastic monotonicity of the contact process stated in ( [23l >, it suffices to prove 
the result for A = S. We consider a modification of where the process resets every 
time an infection travels between C and D, i.e. each time an infection travels between 
C and D, we start again from every site of S infected. The number of resets f for this 
process stochastically dominates T^g. Moreover, using the estimates from Step 2 and 

Step 3, after each reset there is probability at least that the process dies out before 

any infection can travel between C and D. Hence T is dominated by a geometric random 
variable with parameter yielding (|44]). 

Step 5 (Main estimate ([^). We have 

Total number of infections exiting through dS 

E 

for the contact process 

We consider the contact process Let Ti, tj_, , T^c denote the times where an 
infection crosses from C to D and let X\, .. . ,x-pc denote the sites of D where these 

IS 

infections arrive. We can decompose the process inside D as a superposition of 
contact processes i.e. 

?psWnD= U 

Ti<t 


- 2 


1.01 
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Since these contact process are created with the same graphical construction, they are 
not independent however, conditionally on the event | > fc, x, = x|, the process 

— Tj), t > Tj) has the same law as a contact process restricted to D and started 
with only the site x, infected. On the other hand, we can write 


Total number of infections exiting dS 
for the contact process 


where N, is the total number of infections generated by exiting D through dS. Thus, 
combining the results of Step 1 and Step 4, we get 


-1 

rf'c 1 

^ IS 

< E 


- 

k^l 


Total number of infections exiting dS 

:C 


L E E 


fc=l xeD 
d{x,C)=l 


for the contact process ^ ^ 

= E E p{r^fs>t, = 

k=l xeD 
d{x,C)=l 

00 , 

^E E PI > £, 

k=l xeD ^ ^ 

d{x,C)=l 

00 ^ , 




Total number of infections exiting 
dS for the contact process 


k=l 
oo 

^E(i 

k=l 


1 


4r(C) 

< 4r(C)g-1.01r(C)iOi 




“ 2 


1.01 


Step 6 (Main estimate (361). We have 


Extinction time of the contact process 


We use the same idea as in Step 5. First, we write 


< e3Kc)_ 


Extinction time of 


POO 


+ E 


/o l{ePs(OnD^0} 


dt 


Recalling the notation of Sfep 5, and denoting by fj the extinction time of f we have 


poo 


1 


< E 

Ef^ 

k=l 
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Copying the arguments we used in the previous step, we find that 


(3 


k=l xeD 
d{x,C)=l 



< 2c4''(^). 


X 


E 


Extinction time of 


( 45 ) 


It remains to boimd the expected total time during when the contact process has an 
infected site in C. Again, we decompose the process inside C as a superposition 
of contact processes. More rigorously, let 71,72 ,..., 7^c denote the times where an 

IS 

infection crosses from D to C and let 1/1,.. • ,i/-y c denote the sites in C where these 

I® 

infections arrive. We write n C as a superposition of contact processes 

7s(t)nc = ffc(0u U 

li<t 


This yields 


/„ lKfs(t)nc#0)'^‘ 


E[r] + y; E 

k=l 


^{^Ls>n 




where T is the extinction time of and T) is the extinction time of Conditionally 


c 


on the event | > fc|, the process is stochastically dominated by Thus, using 

the estimates of Step 3 and Step 4, we get 

PCO 1 00 f / _ >. 

® i - E[T] + EP{rfi>*:}E(T] 

< (46) 

Finally, putting together (|4^ and (|4^, we find that 


Extinction time of ^ 


< 2c4''(‘^) + 


which completes Step 6 and finishes the proof of the main estimates. 


□ 


5 Questions and possible extensions. 

There are several natural questions left open in the paper - some have already been 
stated. Here we present a few more which, in our opinion, might be interesting to look 
at. 

First, we have seen that the critical parameter pc for Bernoulli CMP on d-dimensional 
lattices is non trivial. Can this result be generalized to other graphs? We believe this to 
be true with minimal assumptions on the graphs and make the following conjecture: 
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Conjecture 5.1. Bernoulli CMP on any graph with bounded degrees has a non-trivial phase 
transition. 


A first step to prove this assertion might be to show the result for trees with bounded 
degrees. By coupling, it suffices to consider regular trees, which seems to be a much 
easier problem than the general case. However, it does not seem straightforward to 
extend the result to general graphs since (contrarily to the contact process for example) 
we cannot directly compare Bernoulli CMP on a given graph to Bernoulli CMP on a 
universal cover of the same graph. 

A similar question applies for Continuum CMP: given a graph G, what are the 
conditions on the distribution of the radii r to ensure the existence of a sub-critical 
phase? If G has exponential growth, it is clear that r should, at least, admit some 
exponential moments. Is this sufficient, at least tor trees? 

Concerning the model of degree-weighted CMP, an important setting is that of 
Gabon-Watson trees for which we expect: 


Conjecture 5.2. Let a >\ and let G be a Galton-Watson tree with reproduction law B such 
that E[exp(cB“)] < oo for any c > 0. Then, degree weighted CMP on G with expansion 
exponent a has a non-trivial phase transition. 


In view of Theorem 4.1 this conjecture, if true, implies that the contact process on 
Galton-Watson trees has a non trivial phase transition whenever the progeny distribu¬ 
tion B has very light tails. This is a work in progress. 


There are also many questions regarding finer percolation properties of the GMP 
which might be interesting to study. For example, is there percolation at criticality? 
When p 7 ^ pc, what is the tail distribution of the size of finite clusters? Is the decay 
faster than exponential? 


Another possible direction of investigation is to consider more general definitions 
of the GMP. As stated at the end of Section]^ it is possible to generalize the notion of 
admissibility. For example, what happens if we consider an expansion exponent smaller 
than 1. Do we have Pc(^) = 1 foi" Bernoulli GMP whenever a; < 1? 


Finally, concerning the connection between the GMP and the contact process, it 
would be extremely satisfying to extend Theorem |4.1| to a = 1. Assuming Gonjecture 

any Galton-Watson tree whose progeny distribution admits exponential moments of all 
orders. This is in particular the case of the Poisson distribution appearing in the limit 
of Erdos-Renyi random graphs. Gonversely, if degree weighted GMP always has an 
infinite cluster for a = 1, does this imply that the critical infection rate for the contact 
process is zero? 

Acknowledgments. A.S. would like to thank J.-B. Gouerefor stimulating discussions con¬ 
cerning continuum percolation and related models. 


5.2 this would imply the existence of a sub-critical phase for the contact process on 
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